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Abstract. We review the recent progress in the theoretical description and experimental observation of 
multiple parton interactions. Subjects covered include experimental measurements of minimum bias in- 
teractions and of the underlying event, models of soft physics implemented in Monte Carlo generators, 
developments in the theoretical description of multiple parton interactions and phenomenological studies 
of double parton scattering. This article stems from contributions presented at the Helmholtz Alliance 
workshop on ” Multi-Parton Interactions at the LHC”, DESY Hamburg, 13-15 September 2010. 


PACS. 11.80.La Multiple scattering — 12.38.Bx Perturbative calculations — 12.38.Lg Other nonpertur- 
bative calculations — 12.38.Qk Experimental tests — 12.39.St Factorization — 13.87.-a Jets in large-Q2 
scattering — 13.87.Fh Fragmentation into hadrons — 14.70.Fm W bosons 
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1 Introduction 


The Large Hadron Collider (LHC) began operation in 2008, opening a new chapter in particle physics. The basis for 
understanding hadronic collisions at high energy is provided by the QCD improved parton model. In this framework 
each hadron is described as a collection of essentially free elementary constituents. The interactions between con- 
stituents belonging to different colliding hadrons are the seeds of the complicated process which eventually leads to 
the particles observed in the detector. Due to the composite nature of hadrons, it is possible to have multiple parton 
hard-scatterings, i.e. events in which two or more distinct hard parton interactions occur simultaneously in a single 
hadron-hadron collision. At fixed final state invariant masses, such cross sections tend to increase with collision energy 
because partons with successively lower momentum fraction x, hence rapidly increasing fluxes, are being probed. As a 
result, events with relatively low invariant masses could receive enhanced contributions from multiple hard scatterings. 
This class of events is known as Multiple Parton Interactions (MPI), while those in which only a single pair of partons 
produce a hard scattering are referred as Single Parton Scattering (SPS). 

The MPI can manifest themselves in various ways in high energy hadronic collisions. It is natural to expect a 
relation between the multiplicity of simultaneous partonic scatterings and their typical scale. In particular, large 
hadronic activity is observed in the soft regime, characterized by small transverse momenta (pr) of the produced 
particles. For relatively large pr values, the observation of MPI will mostly focus on two simultaneous scatterings, i.e. 
on Double Parton Scattering (DPS). Unfortunately, also depending on the scale of the partonic subprocess our ability 
to describe MPI in the pQCD framework is challenging. Whereas it is most legitimate to use pQCD methods for the 
description of MPI at large pr, it is necessary to supplement the pQCD picture in the soft regime with models of soft 
physics . 

The evidence for MPI comes from high pr events observed in hadron collisions at the ISR at CERN [I] and later 
at the Fermilab Tevatron collider [2\[3|4]. At lower pr, underlying event (UE) observables have been measured in pp 
collisions in dijet and Drell-Yan events at CDF in Run I [5] and Run II [6] at center-of-mass energies of ys = 1.8 TeV 
and 1.96 TeV respectively, and in pp collisions at \/s = 900 GeV in a detector-specific study by CMS [7]. 

At small transverse momentum MPI have been shown to be necessary for the successful description of the UE in 
Monte Carlo generators such as PYTHIA [8//9\[10] or HERw1Ic [11][12]. Additionally, MPI are currently invoked to account 
for observations at hadron colliders that would not be explained otherwise: the cross sections of multi-jet production, 
the survival probability of large rapidity gaps in hard diffraction, etc. [13]. The wide range of phenomena in which 
MPI are involved highlights the urgency of a more thorough understanding of these reactions both experimentally and 
from a theoretical point of view. 

The last few years have proven to be a renaissance for research work on MPI. The renewed interest in the field 
follows from the expected abundance of MPI phenomena at the LHC and thus their importance for the full picture 
of hadronic collisions, as well as opportunities provided by the LHC to measure multiple parton hard-scatterings. In 
particular, given the inability to describe the very soft regime with perturbative methods, the relevance of experimental 
measurements by the LHC collaborations of observables containing information on MPI, sensitive to the underlying 
event or minimum bias events, cannot be overstated. 

Ultimately, one would strive for a uniform and coherent description of MPI in both soft and hard regimes. At 
present we are still far away from this goal, as so far essentially separate research efforts focus on specific aspects of 
MPI. In the present article we attempt to bridge this gap by reviewing the current status of the field. The article is 
a result of the Helmholtz Alliance workshop “Multi-Parton Interactions at the LHC” which took place at DESY in 
September 2010. Its main goal was to bring the experimental and theoretical MPI community together, providing a 
forum for a discussion and scientific exchange. 

This article is organized as follows. We begin with a review of the experimental measurements (available at the 
time of the workshop) in Chapter [2} Progress in the implementation of MPI in the Monte Carlo event generators is 
described in Chapter [3] We then focus on the theoretical aspects of MPI in Chapter fand discuss the phenomenology 
of selected DPS processes at the LHC in Chapter [5] 


2 Experimental situation 


A complete description of hadronic activity in high energy collisions requires understanding of the UE, as it constitutes 
the unavoidable background to most observables. From an experimental point of view, the UE gathers all the activity 
accompanying the actual hard scattering one is interested in measuring. In this sense the UE consists of MPI and the 
interactions between consituents of beam remnants, left behind after the scattering partons have been pulled out. 

Even more inclusive measurements probe the so-called Minimum Bias (MB) events. These are events which are 
collected with a relatively non-restrictive trigger which accepts large fraction of events, of both the elastic and inelastic 
nature, inelastic involving the diffractive as well as soft and hard “hard-core” events. 

Since it is impossible to uniquely separate the UE from the hard scattering process on an event-by-event basis, 
the topological structure of the outcome of hadronic collisions is focused on instead. Typically, studies of UE rely 
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Fig. 1. Definition of regions in the azimuthal angle with respect to the leading track. 


on measurements of the properties of charged particle production, while the charged multiplicity and pr spectrum 
are basic MB observables. Observing charged particles allows one to investigate the region of very low pr, crucial 
for exploring soft and semi-hard physics. At higher values of pr, it is possible to directly observe hard MPI in the 
form of double parton scattering (DPS). In this chapter, we review the experimental studies of UE, MB, forward 
measurements, as well as discuss the prospects for observing hard DPS at the LHC. 


2.1 Underlying Event Measurements at ATLAS 
Contributing author: D. Kar (on behalf of the ATLAS Collaboration) 


This section reports on the measurement of UE observables, performed with the ATLAS detectorat the LHC using 
proton-proton collisions at center-of-mass energies of 900 GeV and 7 TeV [14]. 

At the detector level, charged particles are observed as tracks in the inner tracking system. The direction of the 
track with the largest pr in the event — referred to as the “leading” track — is used to define regions of the 7—-¢ plane 
which have different sensitivities to the UE. The axis given by the leading track is well-defined for all events, and is 
highly correlated with the axis of the hard scattering in high-pr events. A single track is used as opposed to a jet or 
the decay products of a massive gauge boson, as it allows significant results to be derived with limited luminosity and 
avoids the systematic measurement complexities of alignment with more complex objects. 

As illustrated in Fig. |1| the azimuthal angular difference between charged tracks and the leading track, |A¢| = 
|Ø — dreading track|, is used to define the following three azimuthal regions [5]: 


— |A¢| < 60°, the “toward region”; 
— 60° < |A¢| < 120°, the “transverse region”; and 
— |A¢| > 120°, the “away region”. 


The transverse regions are most sensitive to the underlying event, since they are generally perpendicular to the axis 
of hardest scattering. 
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2.1.1 Analysis Details 


All data used in this paper were taken during the LHC running periods with stable beams and defined beam-spot 
values, between 6th—15th December 2009 for the analysis at \/s = 900 GeV, and from 30th March to 27th April 2010 
for the 7 TeV analysis. Only events with leading track pr > 1 GeV within the inner detector, |7| < 2.5, were considered, 
in order to reject events where the leading track selection can potentially introduce large systematic effects, and also 
to reduce the contribution from diffractive hard scattering processes. All the other tracks were required to have pr 
> 500 MeV and the same 7) range. 

The 900 GeV and 7 TeV data respectively correspond to integrated luminosities of 7 pbt and 168 pb™ t, respectively, 
and the effects of pileup was negligible. For the MC models considered here, the contribution of diffractive events to 
the underlying event observables was less than 1%. 

The data were corrected back to charged primary particle spectra satisfying the event-level requirement of at least 
one primary charged particle within pr > 1 GeV and |n| < 2.5. A two step correction process was used, where first the 
event and track efficiency corrections were applied, then an additional bin-by-bin unfolding was performed to account 
for possible bin migrations and any remaining detector effects. 


2.1.2 Results 


In this section, corrected distributions of underlying event observables are compared to model predictions tuned to a 
wide range of measurements. The transverse, toward and away regions each have an area of Ad An = 107/3 in n-¢ 
space, so the density of particles (d? Nen /dndọ) and transverse momentum sum (d?)“pr/dnd¢) are constructed by 
dividing the mean values by the corresponding area. The leading charged particle is included in the toward region 
distributions, unless otherwise stated. 

The data, corrected back to particle level in the transverse, toward and away regions are compared with predictions 
by PYTHIA with the ATLAS MC09 [16], DW [I7], and Perugia0 [18] tunes, by HERWIG+JIMMY with the 
ATLAS MCO09 tune, and by PHOJET [20]. The ratios of the MC predictions to the data are shown at the bottom of 
these plots. The error bars show the statistical uncertainty while the shaded area shows the combined statistical and 
systematic uncertainties. 

The charged particle multiplicity density, is shown in Fig. 

The average number of charged particles in the transverse region increases with leading pr, until it reaches an 
approximately constant “plateau”. All the pre-LHC MC tunes considered show at least 10-15% lower activity than the 
data in the transverse region plateau. The PYTHIA DW tune is seen to be the closest model to data for the transverse 
region. The toward and away regions are dominated by jet-like activity, yielding gradually rising number densities. 
The 900 GeV and 7 TeV ATLAS data show the same trend. The underlying event activity is seen to increase by a 
factor of approximately two between the 900 GeV and 7 TeV data. This is roughly consistent with the rate of increase 
predicted by MC models tuned to Tevatron data. 

The charged particle scalar pr sum, X` pr density, is shown in Fig. |3| The summed charged particle pr in the 
plateau characterizes the mean contribution of the underlying event to jet energies. Again, we can see that pre-LHC 
tunes model CDF data better than ATLAS data. The higher number density implies a higher pr density as well. In 
the toward and away regions, jet-like rising profiles are observed, in contrast to the plateau in the transverse region. 
The toward region includes the leading charged particle, and has a higher X` pr than the away region as there is 
higher probability of high-py particles being produced in association with the leading pr charged particle. In the 
toward region the highest fraction of energy has been allocated to a single charged particle. This implicitly reduces 
the number of additional charged particles in that region, since there is less remaining energy to be partitioned. As a 
result the multiplicity of charged particles is slightly lower in the toward region by comparison to the away region for 
high p'e*¢. The increase of the pr densities in the toward and away regions indicates the extent of the variation in the 
charged fraction of the total energy in each region. 

In Fig.|4} the standard deviation of the charged particle multiplicity and charged particle scalar X` pr densities, are 
shown. The mean and standard deviation of the pr density in the transverse region characterize a range of additional 
energy that jets might acquire if the underlying event were uniformly distributed. The confirmation that the magnitude 
of the standard deviations of the distributions are comparable to the magnitudes of the mean values indicates that a 
subtraction of the underlying event from jets should be done on an event by event basis, rather than by the subtraction 
of an invariant average value. 

The correlation between the mean pr of charged particles and the charged particle multiplicity in that region is 
sensitive to the amount of hard (perturbative QCD) versus soft (non-perturbative QCD) processes contributing to the 
underlying event. 

In Fig. |5| the ATLAS profiles in the transverse and away regions are very similar, showing a monotonic increase 
of (pr) with Nen. The profile of the toward region is different, as it is essentially determined by the requirement of 
a track with pp > 1 GeV. For Nen = 1, it contains only the leading charged particle and as the Nen is increased by 
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Fig. 2. ATLAS data at \/s = 900 GeV (left) and 7 TeV (right), showing the density of the charged particles in the transverse 
region (top row), and in all three regions (bottom row) 


inclusion of soft charged particles the average is reduced. However, for Nen > 5 jet-like structure begins to form, and 
the weak rise of the mean pr is observed. Comparing the 900 GeV and 7 TeV data, it is seen that the mean charged 
particle pr vs. Nen profiles are largely independent of the energy scale of the collisions. 

The angular distributions with respect to the leading charged particle of the charged particle number and X` pr 
densities at the center-of-mass energy of 7 TeV at ATLAS, are plotted in Fig [6] The leading charged particle taken to 
be at A¢ = 0 has been excluded from the distributions. The data are shown for four different lower cut values in leading 
charged particle pr. These distributions are constructed by reflecting |A@| about zero, i.e. the region —r < Ad < 0 is 
an exact mirror image of the measured |A¢| region shown in 0 < Ad < z. 

These distributions show a significant difference in shape between data and MC predictions. With the increase of 
the leading charged particle pr, the development of jet-like structure can be observed, and the corresponding sharper 
rise in transverse regions compared to the MC. MC models essentially predict a stronger correlation than is seen in 
the data, and this discrepancy in toward region associated particle density was also observed at CDF [21]. 

A complementary way [22] to look at the angular correlation is by either subtracting the minimum of the distribution 
(determined by a second-order polynomial fit), or by subtracting the opposite side distribution (defined according to if 
pseudorapidity has the same or the opposite sign as the leading track) from the same side distribution and normalizing 
to unity. In Fig.|7| it is seen that the models are better at lower 7 than at higher. 


2.1.3 Summary and Conclusions 


One of the goals of these analyses is to provide data that can be used to test and improve MC models for current and 
future physics studies at the LHC. The underlying event observables presented here are particularly important for 
constraining the energy evolution of multiple partonic interaction models, since the plateau heights of the underlying 
event profiles are highly correlated to multiple parton interaction activity. The data at 7 TeV are crucial for MC 
tuning, since measurements are needed with at least two energies to constrain the energy evolution of MPI activity. 
There is no current standard MC tune which adequately describes all the early ATLAS data. However, using 
diffraction-limited minimum bias distributions and the plateau of the underlying event distributions presented here, 
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Fig. 3. ATLAS data at \/s = 900 GeV (left) and 7 TeV (right), showing the scalar X` pr density of the charged particles in 
the transverse region (top row), and in all three regions (bottom row) 


ATLAS has developed a new PYTHIA tune AMBT1 (ATLAS Minimum Bias Tune 1) and a new HERWIG+ JIMMY tune 
AUET1 (ATLAS Underlying Event Tune 1) which model the pr and charged multiplicity spectra significantly better 
than the pre-LHC tunes of those generators [2324]. It is critical to have sensible underlying event models containing 
our best physical knowledge and intuition, tuned to all relevant available data. 


2.2 Multiple Parton Interactions Studies at CMS 
Contributing authors: P. Bartalini and Livio Fand (on behalf of the CMS Collaboration) 


This section summarizes the early Underlying Event and forward measurements of the CMS collaboration using pp 
collision data up to highest energies of ys = 7 TeV. It also reports along the feasibility study for the direct measurement 
of double parton scattering phenomena focusing on the 3jet + y channel. 

A detailed description of the CMS detector is available in Ref. [25]. Generator level Monte Carlo (MC) predictions 
are compared to the data corrected with a bayesian unfolding technique taking into account the detector effects ; 

The predictions for inelastic events are provided here by several tunes of the PYTHIA program, versions 6.420 [27] 
[15] and 8.134] [28129]. PHOJET [30] is also used in the forward measurements: 

The pre-LHC tune D6T [BI13) of Pyra 6, which describes the lower energy UA5 and Tevatron data, is a widely 
used reference that will also be used for most of the presented analyses. The tunes DW [13] and CW [7], which were 
found to describe best the UE CMS data at 0.9 TeV whereas D6T predictions were too low [7], will also be discussed 
for the 7 TeV data. The pre-LHC tune Perugia-0 and the new tune, Z1 [33], adopt pr ordering of parton showers 
and the new PyTHIA MPI model [34]. It includes the results of the Professor tunes [35] considering LEP fragmentation 
and the color reconnection parameters of the AMBT1 tune [36], while with the first CMS UE results [737] have been 
used to tune the parameters governing the value and the ys dependence of the cut-off transverse momentum that 
in PYTHIA regularizes the divergence of the leading order scattering amplitude as the final state parton transverse 
momentum pr approaches 0. The tune Z2 is similar to Z1, except for the transverse momentum cut-off at the nominal 
energy of \/s9 = 1.8 TeV which is decreased from 1.932 GeV/c to 1.832 GeV/c. PYTHIA 8 also uses the new PYTHIA 


1 PYTHIA version 8.108 is used in the feasibility studies reported in section 
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Fig. 4. TLAS data at ys = 900 GeV (left) and 7 TeV (right), showing the standard deviation of the density of the charged 
particles (top row) and the standard deviation of the scalar }> pr density of charged particles (bottom row) in the transverse 
region. 
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Fig. 6. ATLAS data showing the ¢ distribution of charged particle multiplicity (left) and scalar X` pr density (right), with 
respect to the leading charged particle rotated to ¢icading = 0, excluding the leading charged particle and compared to different 
MC model predictions. The distributions obtained by restricting the minimum leading charged particle pr to different values 
are shown separately. The plots were symmetrized by reflecting them about Ad = 0 
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Fig. 7. The Ag crest shape obtained by subtracting the minimum are shown in (a) and (c), while those obtained by the 
subtracting the ‘opposite from same’ are shown in (b) and (d). The left two plots are for |n| < 1.0 and the right two plots are 
for |n| < 2.5. 


MPI model, which is interleaved with parton showering. The default Tune 1 is adopted here. PYTHIA 8 includes 
soft and hard diffraction [38], whereas only soft diffraction is included in PYTHIA 6; the diffractive contributions 
are, however, heavily suppressed by the trigger and event selection requirements, especially for large pr values of the 
leading track-jet. The parton distribution functions used to describe the interacting protons are the CTEQ6LL set for 
D6T and Z2 and the CTEQSL set for the other simulations [B940]. 


2.2.1 The Early Underlying Event Measurements 


In the presence of a hard process, characterized by particles or clusters of particles with a large transverse momentum pr 
with respect to the beam direction, the final state of hadron-hadron interactions can be described as the superposition 
of several contributions: products of the partonic hard scattering with the highest pr, including initial and final state 
radiation; hadrons produced in additional MPI; “beam-beam remnants” (BBR) resulting from the hadronization of 
the partonic constituents that did not participate in other scatterings. Products of MPI and BBR form the UE, which 
cannot be separated from initial and final state radiation. 

The early CMS UE measurements focus on the understanding of the UE dynamics studying charged particle 
production with two different approaches. The first (traditional) approach [7[37] concentrates on the study of the 
transverse region, which is defined considering the azimuthal distance of the reconstructed tracks with respect to the 
leading track or leading track-jet of the event: 60° < |Ad| < 120°. The jet reconstruction algorithm used in these 
studies is SisCone [41]. On top of the traditional approach, a new methodology using anti-kr jets and relying 
on the measurement of their area [43] is adopted for the first time by CMS using charged particles in pp collision 
data collected at ys = 0.9 TeV [44]. The new set of UE observables consider the whole pseudorapidity-azimuth plane 
instead of the transverse region and inherently take into account the leading jets of an event. 
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Fig. 8. (Upper plots) average multiplicity and average scalar ` pr in the transverse region as a function of the leading track- 
jet pr, for data at ys = 0.9 TeV and ys = 7 TeV. (Bottom left plot) normalized scalar X` pr distribution in the transverse 
region for data at ys = 0.9 TeV and ys = 7 TeV; the leading track-jet is required to have pr > 3 GeV/c. Predictions from 
PYTHIA 6 tune Z1 and PYTHIA 8.135 Tune 1 are compared to the corrected data. The inner error bars indicate the statistical 
uncertainties affecting the measurements, the outer error bars thus represent the statistical uncertainties on the measurements 
and the systematic uncertainties affecting the MC predictions added in quadrature. (Bottom right plot) normalized median 
of pr over jet area for track-jets reconstructed from collision data at ys = 0.9 TeV (black circles). Predictions from several 
PYTHIA 6 tunes and PYTHIA 8 Tune 1 are compared to data. 


The centre-of-mass energy dependence of the hadronic activity in the transverse region is presented on the two 
top Figures [8] as a function of the pr of the leading track-jet. The data points represent the average multiplicity and 
average scalar track-pp sum dependence, for ys = 0.9 TeV and ys = 7 TeV using tracks with a pseudorapidity 
|n| < 2.0 and pr > 0.5 GeV/c. A significant growth of the average multiplicity and of the average scalar pr sum of 
charged particles transverse to that of the leading track-jet is observed with increasing scale provided by the leading 
track-jet pr, followed by saturation at large values of the scale (more evident for multiplicity profile than average 
scalar pr sum). A significant growth of the activity in the transverse region is also observed, for the same value of the 
leading track-jet pr, from ys = 0.9 TeV to ys = 7 TeV. These observations are consistent with the ones obtained at 
Tevatron [45]. The evolution with the hard scale of the ratio of the UE activity at 7 TeV and 0.9 TeV is remarkably well 
described by the Z1 MC. The trend is also reproduced by PYTHIA 8. The Z2 predictions at ys = 0.9 TeV (not shown 
here) agree with Z1 in shape, but the normalization is 5-10% higher for both the observables; this trend is opposite 
with respect to the one observed at 7 TeV and indicates that a less pronounced ,/s dependence of the transverse 
momentum cut-off should be adopted for tunes using the CTEQ6LL PDF set than for the tunes optimized for the 
CTEQSL set. 
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The strong growth of UE activity with charged particles is also striking in the comparison of the normalized 
distributions of charged particle multiplicity (not shown here) and of scalar pr sum which is presented in bottom-left 
plot of Figure [8] for events at /s = 0.9 TeV and ys = 7 TeV with leading track-jet pr > 3 GeV/c. The particle 
pr spectrum (not shown) extends up to pr > 10 GeV/c, indicating the presence of a hard component in particle 
production in the transverse region. The distributions for track-jet pr > 3 GeV/c, which extend up to quite large 
values of the selected observables in the transverse region are quite well described by the various MC models, over 
several orders of magnitude. This observation gives support to the implementation of MPI in PYTHIA. 

The novel technique to quote the UE activity relies on the introduction of “ghosts”, virtual deposits of very low 
energy filling the overall phase space which are taken into account by the jet clustering algorithm. The estimator of 
the overall soft background activity in an event can be derived as the median of the ratio between the transverse 
momentum and the area of the jets. One of the advantages of using the median compared to the mean is that it turns 
out to be less sensitive to the influence of outliers, i.e. in particular the leading jets in an event. In order to cope with 
the low occupancy observed at ys = 0.9 TeV, CMS redefines such observable restricting the median only to those jets 
which have physical deposits on top of ghosts: 


p! = median eal -C (1) 


where C is the occupancy of the event, which is the summed area )> j A; covered these jets divided by the considered 
detector region Atot. In the CMS analysis at ys = 0.9 TeV, jets are reconstructed with the anti-Kr algorithm [42] 
using tracks with |ņ| < 2.0 and pr > 0.3 GeV/c. In the right bottom plot of Figure [B] the p' observable is presented 
for minimum bias events. The general pattern of deviations from data with respect to the considered PYTHIA tunes 
looks rather similar to the one observed with the traditional UE measurement. 


2.2.2 Study of the Activity in the Forward Region 


CMS reports a measurement of the energy flow in the forward region (3.15 < |n| < 4.9, where 7 denotes the pseudora- 
pidity) for minimum bias and dijet events in pp interactions with centre-of-mass energies \/s of 0.9 TeV, 2.36 TeV 
and 7 TeV. This measurement is connected to the ones reported in the previous sections as the basic philosophy is the 
same: it concentrates on the complementary activity of a pp interaction for different energy scales of the reconstructed 
leading objects. 

The energy flow in the region of the Hadron Forward detector is measured in two different event classes: in minimum 
bias events and in events with a hard scale provided by a dijet system at central pseudorapidities (|ņ| < 2.5) and with 
transverse energy Er, jet > 8 GeV for ys = 0.9 TeV and 2.36 TeV; the dijet threshold is increased to 20 GeV for ys = 
7 TeV. The results are qualitatively similar at all the investigated centre of mass energies. Fig. [9] shows the results of 
the forward energy flow at \/s = 7 TeV for the two event classes compared to predictions from Monte Carlo event 
generators. The measured forward energy flow is found to be significantly different between the two event classes, with 
a sensitive increase and a more central activity seen in dijet events. 


2.2.3 Multiple Parton Interactions in High-pr Phenomenology 


Quantifying the MPI cross sections basically deals with the measurement of oe ff, the scale factor which characterizes 
the inclusive rate of the interactions [47/48], cf. Chapters |4| and |5| From a phenomenological point of view oes f is 
a non perturbative quantity related to the transverse size of the hadrons and has the dimensions of a cross section. 
The measurements performed by the AFS, CDF and DO collaborations [B950] favor smaller values of o. ¢¢ with 
respect to the naive expectations. The consequently increased rates of multiple parton interactions can be interpreted 
as an effect of the hadron structure in transverse plane [51]. Extending such measurements at the LHC and studying 
the possible scale dependency of efr is definitely of great interest and may have a deep impact on the data driven 
estimations of the MPI backgrounds to new physics. 

The production of four high-pr jets is the most prominent process to search for multiple high pr scatterings: two 
independent scatters (i.e. DPS) in the same pp or pp collision, each producing two jets. Such a signature has been 
searched for by the AFS experiment at the CERN ISR, by the UA2 experiment at the CERN SppS and by the CDF 
and DO experiment at the Fermilab Tevatron. 

However, searches for double-parton scattering in four-jet events at hadron colliders may face significant back- 
grounds from other sources of jet production, in particular from QCD bremsstrahlung (Fig.|10}left). Typical thresholds 
employed in jet triggers bias the event sample towards hard scatterings. However, a high-pr jet parton is more likely 
to radiate additional partons, thus producing further jets. Thus, the relative fraction of jets from final-state showers 
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Fig. 9. Energy flow in the minimum bias (left) and di-jet (right) samples as a function of Nobs at \/s = 7 TeV. Uncorrected data 
are shown as points, the histograms correspond to the MC predictions. Error bars corresponds to statistical errors. The shaded 
bands in these plots represent the systematic uncertainties of the measurements, which are largely correlated point-to-point. 
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Fig. 10. Definition of azimuthal angle between pairs, together with typical configurations of double-bremsstrahlung (left) and 
double-parton scattering events (right). 


above a given threshold is enlarged in jet trigger streams which is an unwanted bias. On the other hand, looking for 
four jets in a minimum-bias stream will yield little statistics. 

Therefore the strategy to directly measure the MPI rate in high-pr regime at hadron colliders also includes the 
study of multi-jets or jets+photon final states. Indeed the CDF and DO collaborations studied final states with one 
photon and three jets looking for pairwise balanced photon-jet and dijet combinations. The data sample was selected 
with the experiment’s inclusive photon trigger, thereby avoiding a bias on the jet energy. The better energy resolution 
of photons compared to jets purifies the identification of Er balanced pairs. Tevatron found an excess in pairs that 
are uncorrelated in azimuth with respect to the predictions from models without additional hard parton scatters per 
proton-proton scatter. CDF interpreted this result as an observation of double-parton-scatters. 

Analyses trying to identify two hard scatters in multi-jet events typically rely on methodologies which overcome 
combinatorics. There are three possible ways to group four objects into two pairs: combinations are commonly selected 
pairing objects which are balanced in azimuth and energy. The flavor or other specific features of the jets may be used 
to decrease the combinatorics and to make looser the constraints on the balancing. One example of such a final state 
is constituted by events with two b jets and two additional light jets. 

In order to discriminate double-parton scatters against double-bremsstrahlung events, CMS studies the observables 
Ao, employed by AFS, and Ad), employed by CDF, probing the azimuthal angle between pairs (Fig. [1o}. 
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Fig. 11. Differential cross section shape as a function of Ad (upper plots) and Ad (bottom plots) variables. Predictions 
from PYTHIA 8.108 (Default scenario) and HERWIG 2.2.0 (left panel) and from three different PYTHIA settings (right panel) 
shown. 


Expectations for the above described variables are therefore Ao?) = m /2 and Ad) = x if additional jets come from 
double-bremsstrahlung. Otherwise, i. e. if additional jets come from multiple interactions, both variables should be 
distributed uniformly. 

Differential cross section shape predictions for the A@ observables in pp interactions at 14 TeV are shown in 
Fig. HERwIG 2.2.0 [19/52] and PYTHIA 8.108 with default settings which include multiple interactions and 
showering predict similar shapes (Fig. |11}left). The discrimination power of the selected observables to Multiple 
Parton Interaction patterns is clearly shown in Fig. right, where events with MPI switched off (Shower scenario) 
are compared to events with parton shower switched off (MI scenario). The differences are particularly pronounced 


when selecting the Ad‘) observable. 
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2.2.4 Conclusions 


A strong growth of the UE activity is observed with increasing leading track-jet pr for both ys = 7 TeV and ys = 0.9 
TeV. At 7 TeV this fast rise is followed above ~8 GeV/c by a saturation region with nearly constant multiplicity 
and small pr increase. The same pattern is observed at 0.9 TeV, with the saturation region starting at ~4 GeV/c. A 
strong growth of the activity is also observed with increasing centre-of-mass energy. The large increase of activity in 
the transverse region is also observed in the ` pr distribution, indicating the presence of a hard component in the 
transverse region. Very good post-LHC MC tunes are available for the description of the UE in the central region. 

A measurement of the underlying event using the jet-area/median approach is also reported, demonstrating its 
sensitivity to different underlying event scenarios. 

Complementary underlying event measurements in the forward region are also presented. The energy flow in the 
forward direction is measured for minimum bias and central di-jets events. A more global UE description including 
both the central and the forward regions is certainly one of the next MC tuning challenges, with deep impact on the 
understanding of the MPI dynamics. 

The Multiple Parton Interactions measurement strategy in the high-pr regime is also briefly discussed focusing 
on the 3jet + y topology. The very good performances of the LHC machine should allow to have soon the integrated 
luminosity conditions adequate to perform these measurements over a wide range of energy scales, with deep impact 
on the data driven estimation of the MPI backgrounds to searches. 


2.3 Minimum Bias Physics at LHCb 
Contributing author: M. Schmelling (on behalf of the LHCb Collaboration) 


At the startup in 2009 the LHC provided proton-proton collisions at a center-of-mass energy of ys = 0.9TeV. 
Although higher collision energies have previously been reached at proton-antiproton colliders, it was for the first 
time that the TeV-scale was studied in proton-proton collisions. Using a data sample with an integrated luminosity 
of only 6.8 + 1.0 ub~! recorded by the LHCb detector, a first measurement of the production cross-section of neutral 
K-mesons was performed in a kinematic range not accessible to the other LHC experiments. In 2010 the collision 
energy was moved up to 7 TeV and the performance of the machine improved exponentially. Until September LHCb 
collected an integrated luminosity of over 3 pb~! with a data acquisition efficiency larger than 90%. At the same time 
the detector calibration approached its design values. In the following, after a brief description of the detector the 
first measurements on strangeness production and studies of baryon number transport and baryon suppression in the 
fragmentation will be presented, before finally discussing prospects for doing diffractive physics with LHCb. 


2.3.1 The LHCb Experiment 


The LHCb detector [53] is a forward spectrometer, covering the angular range of 15 < 0 < 300mrad with respect 
to the beam axis. A schematic view of the experiment is shown in Fig [12] The detector offers tracking, calorimetry 
and particle identification over most of its forward acceptance. Momenta of charged particles are determined from 
the deflection by a dipole magnet with a field integral of 4Tm. The interaction region is surrounded by the Vertex 
Locator (VeLo). Going downstream, a first RICH detector and the so-called TT tracking station are still located in 
front of the magnet. Immediately behind the magnet follows the second part of the tracking system, consisting of 
a high granularity Inner Tracker (IT) in the region of large particle densities close to the beam pipe and the Outer 
Tracker system at larger transverse distances. VeLo, TT and IT are silicon strip detectors, the OT consists of straw 
tubes. Following the tracking system is a second RICH detector, a pre-shower and scintillating pad detector (SPD/PS), 
electromagnetic calorimeter (ECAL), hadron calorimeter(HCAL) and muon system for the identification of electrons 
and photons, neutral hadrons and muons, respectively. The RICH detectors allow pion, kaon, proton separation in the 
momentum range between 2 < p < 100 GeV/c. 

The VeLo has 21 double-layer sensor planes for measuring space points and two single-layer planes providing only 
radial track coordinates. Its layout is shown in Fig.}13| The angular acceptance is larger than for the rest of the 
tracking system and covers also part of the backward hemisphere. However, being located outside of the magnetic 
field, VeLo track segments do not have momentum information. Furthermore, since at least three planes are required 
to reconstruct a track segment, the VeLo is blind in the central region. Charged particle tracks are reconstructed in 
the rapidity ranges —4 < ņ < —1.5and15<7 < 5. 

With the data recorded until summer 2010 a first precise calibration of the various subdetector systems was 
performed. With the tracking system, for example, impact parameter resolutions for high-pr particles around 16 wm 
were reached. Tracking efficiencies above 95% for charged particles with transverse momenta above pr = 200 MeV/c 
were determined by a tag and probe approach using daughter particles from K9 and J/w decays. The measurements 
were found to be in good agreement with Monte Carlo simulations, allowing to set a limit on the systematic uncertainty 
of the tracking efficiency to 3% for high momentum tracks and 4% for soft particles. 
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Fig. 12. Schematic view of the LHCb single arm forward spectrometer. The interaction region is located on the left inside the 
Vertex Locator. The tracking system and the RICH detectors for particle identification are installed both before and after the 
dipole magnet, calorimetry and the muon system are located downstream of the magnet. 
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Fig. 13. Layout of the LHCb Vertex Locator (VeLo) in the horizontal (x, z)-plane, with the z-axis along the direction of the 
proton beams. 21 sensor planes measure space points, the two most backward (—z) layers provide only radial coordinates of 
charged particle tracks. 


2.3.2 Measurements of Strangeness Production 


The study of strangeness production is of particular interest for early measurements since, compared to heavy flavours, 
the cross section is large, and V°-decays (K°, A, A) have a very clean experimental signature which allows to identify 
them unambiguously using only kinematical information even when only a coarse calibration of the detector exists. 
Furthermore, since strange quarks are not present as valence quarks in the initial state and have a mass in an 
intermediate range where QCD predictions have large uncertainties, they directly probe the mechanism of multi- 
particle production in high energy collisions where the theory is least well understood. 


Results from the first cross-section measurements of prompt K°-production in proton-proton collisions at a center- 
of-mass energy ys = 0.9 TeV are shown in Fig.|14} Uncorrelated errors, mainly due to finite statistics and the 
modeling of the shape of the differential cross section in the determination of the correction factors, range between 
6 and 28%. The correlated errors are between 16 and 23% and are dominated by uncertainties in the track finding 
efficiencies and the luminosity measurement. Comparing experimental results and Monte Carlo predictions, one clearly 
sees that the data favour a harder pr-spectrum than the models that were considered. The best description is obtained 
with the Perugia0-tune of PYTHIA, which does not include diffractive processes. 
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Fig. 14. Differential cross-section for K2-production in pp collisions at ys = 0.9 TeV as a function of transverse momentum pr 
and rapidity y. The vertical bars on the data points show the total uncertainties, the purely statistical errors are indicated by 
the tick marks. The histograms are predictions from different settings of the PYTHIA model. The lower plots show the MC/data 
ratios, with the shaded bands the experimental uncertainty. 
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Fig. 15. A/A cross-section ratio as a function of rapidity for center-of-mass energies \/s = 0.9 TeV and \/s = 7 TeV compared 
with expectations from Monte Carlo models (left), and as a function of the rapidity distance to the beam particle (right). 


2.3.3 Baryon Number Transport and Baryon Suppression 


Other observables probing the dynamics of particle production in high energy hadron collisions are cross-section ratios, 
where luminosity and many systematic uncertainties cancel. Results from the study of the A/A cross-section ratio for 
center-of-mass energies \/s = 0.9 TeV and ys = 7 TeV are shown in Fig/15| While A-baryons have two of their three 
valence quarks in common with the proton, all three antiquarks of the ave to be produced in the collision. The 
ratio of the production cross-sections thus measures the baryon-number transport from the beam particles to the 
final state. In general one observes that the measured ratio is lower, i.e. larger baryon number transport, than the 
expectation from the Monte Carlo models. The effect becomes stronger in the forward region and when going to lower 
center-of-mass energies. 

A natural variable to study baryon number transport is the rapidity difference to the beam. While for \/s = 0.9 TeV 
the beam rapidity is at y = 6.6 it rises to yẹ = 8.3 at ys = 7TeV. At the lower energy the LHCb acceptance thus 
reaches much closer to the beam than at \/s = 7 TeV. Plotting the A/A cross-section ratio, as a function of Ay = y,—y, 
also shown in Fig.[15] scaling is observed. The LHCb results are in good agreement with measurements by the STAR 
collaboration. 
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Fig. 16. A/K®& cross-section ratio as a function of rapidity for center-of-mass energies ys = 0.9 TeV and ys = 7TeV. The 
measurements are compared with expectations from Monte Carlo models. (Note the slightly different vertical scales in the two 
plots.) 
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Fig. 17. Schematic classification of inelastic proton-proton collisions. The upper row shows some born-level type diagrams for 
different classes of interactions, the lower row illustrates the angular range into which particles produced in the collision are 
emitted. Note that in the upper row rapidity runs from top to bottom, while it goes from left to right in the lower row. 


Another ratio to look at is the A/K cross-section ratio. Since the A has no valence quarks in common with the 
initial state protons, this ratio measures the suppression of baryon- relative to meson-production. Experimental results 
compared to Monte Carlo predictions are shown in F ig [16] Both data and Monte Carlo show a slight increase in the 
ratio when going from ys = 0.9TeV to ys = 7TeV. This is plausible, since particle masses and kinematic factors 
in general should become less important at higher energies. It is, however, striking that in both cases the ratio is 
significantly underestimated by the Monte Carlo models. Since both the A and the K$ contain a single (anti)strange 
valence quark, the discrepancy cannot be explained by a mismatch in strangeness suppression between data and Monte 
Carlo. Instead it probes directly the understanding baryon formation in the fragmentation process. 


2.3.4 Prospects for the Study of Diffractive Interactions 


A schematic view of different types of inelastic pp-collisions contributing to minimum bias interactions is presented in 
Fig.{17] Here the basic distinction is colour-octet (gluon) and colour-singlet (pomeron) exchange between the colliding 
protons. In the language of QCD the pomeron is understood as a color-singlet two-gluon state. Colour exchange 
implies that the structure of both protons is resolved with the consequence that the colour fields stretched between the 
partons lead to particle production in the full rapidity range. In contrast, colour-singlet exchange is phenomenologically 
described by pomerons coupling to the protons as a whole. No colour is transferred and the protons can either scatter 
elastically or be excited into a high mass state which then decays to produce multi-particle final states. Depending on 
whether only one or both protons are excited these processes are referred to as single- or double-diffractive scattering. 
An example for a higher order process involving pomerons is double pomeron exchange, where both protons stay intact 
and the two pomerons interact to form a massive central system. 
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Fig. 18. Generator level MC study: Mix of single diffractive, double diffractive and non-diffractive events generated by 
PYTHIA 8.135. Here SD1 refers to single diffractive scattering where the excited proton travels in the direction of the LHCb 
detector, in SD2 the decaying heavy mass moves into the opposite direction. The left hand plot shows the mix of events generated 
(black), passing the VeLo track-segment trigger(red) and the diffraction selection (blue). The right hand plot shows the selection 
efficiencies for the two cases. 
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Fig. 19. Generator level MC study: comparison of generated and observable transverse momentum spectra of charged particles 
in the pseudo rapidity range 2 < 7 < 5 when asking for a single VeLo track (left) or, in addition, a rapidity gap (right). 
The spectra are normalized to the number of accepted events. The losses at small transverse momenta are due to incomplete 
geometric coverage at low pr and low 7. 


The above classification of pp-interactions is most adequate at small momentum transfers. With increasing mo- 
mentum transfer, however, the simple picture breaks down. The two protons are resolved into an increasing number of 
partons and the interaction is described by ladder-diagrams of all possible topologies. Diffractive and non-diffractive 
scattering is no longer an unambiguous classification, and even the notion of colour singlet exchange becomes frame 


dependent 

Another conceptual problem with Fig [I 7]is that the upper row represents Born-level type diagrams, i.e. amplitudes, 
while the lower one is a pictorial representation of the cross-section. In many Monte Carlo models, such as e.g. PYTHIA 
[5], the different contributions to the total cross-section are generated independently and interference terms are 
ignored. Furthermore, parameter tuning generally allows to trade e.g. a larger diffractive component against a smaller 
non-diffractive part. It follows that event classification into diffractive and non-diffractive parts based on a Monte 
Carlo implementation is to some extent arbitrary and always model dependent. A better approach would be to avoid 
any such artificial classifications, and instead perform measurements subject to experimental cuts which enhance or 
suppress certain contributions to the cross section. 


To study the prospects for measuring the properties of events with dominantly diffractive contributions, a simple 
generator level analysis, ignoring finite resolution and imperfect efficiencies, of charged particle transverse momentum 
spectra has been performed. The study is based on PYTHIA 8.135 available from [55]. Single proton-proton collisions 
with a center-of-mass energy ys = 7 TeV were generated with process selection pythia. readString("SoftQCD:all=on"). 
The VeLo was simulated with its nominal geometry, and a track was assumed to be measured if at least three stations 
were hit. 
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The event selection was based on the VeLo track segments only. Denoting by ng and nr the number of VeLo 
track segments in the backward (7 < 0) and forward (7 > 0) hemispheres, two selection criteria were studied: (a) 
np +np > 0 and (b) ng == 0 && np > 0. Both criteria ask for at least one track segment in the VeLo, the second 
corresponds to the additional requirement of a rapidity gap An = 2.5 for charged tracks and thus is expected to 
enhance the fraction of diffractive events in the sample. To obtain a quantitative measure for the level of enrichment 
which can be achieved, the PYTHIA process type was analyzed for all events. While evidently giving model dependent 
estimates for the fractions of different events, the qualitative picture is expected to be generic. 

For accepted events the transverse momentum spectrum then was determined using all tracks with a VeLo-segment 
and within the acceptance of the tracking system behind the magnet. The latter was approximated by the requirement 
p > 2GeV/c the pseudo-rapidity range 2 < 7 < 5. Results are shown in Figs [18] and The left hand plot of 
Fig [18] shows the mix of single diffractive, double diffractive and non-diffractive events generated by PYTHIA 8.135. 
One clearly sees that the requirement of a rapidity gap in the backwards region almost completely suppresses non- 
diffractive events while keeping between 20% and 30% of diffractive interactions. The comparison of the transverse 
momentum spectra in Fig [19] shows very good agreement between generated and observed distribution, i.e. a robust 
measurement comparing the fully inclusive transverse momentum spectra and the spectra in events dominated by 
diffraction seems feasible. Measurements of production cross-sections for identified particles and particle ratios are a 
natural extension of these studies. 


2.3.5 Summary and Outlook 


First results from the study of minimum bias events by the LHCb experiment at center-of-mass energies \/s = 0.9 TeV 
and ys = 7 TeV have been presented. Production cross-sections for Kg-mesons at ys = 0.9 TeV were found to have 
harder transverse momentum spectra than expected from Monte Carlo models. The baryon number transport from 
the beam particles to the final state was found to scale with rapidity difference to the beam particles and to be more 
pronounced that expected from the currently used models. The models also feature a stronger baryon suppression 
in the fragmentation than is observed in the experiment. Finally, a model independent approach towards the study 
of diffractive particle production in minimum bias events has been presented. Generator level Monte Carlo studies 
suggest that asking for a rapidity gap of An = 2.5 in the backward region of the VeLO allows to select event samples 
dominated by diffractive processes which, making use of the excellent tracking and particle identification capabilities 
of the LHCb detector, then can be studied in detail in the pseudo-rapidity range 2 < 7 < 5. 


3 Multi-parton Interactions in Event Generators 


The description of low pr hadronic activity used in experimental analyses relies on models implemented in Monte 
Carlo (MC) event generators. These generators combine a perturbative description in terms of multiple scatterings 
with phenomenological models for soft, non-perturbative physics. Recently, the implementation of MPI effects in Monte 
Carlo models has quickly progressed through an increasing level of sophistication and complexity that has deep general 
implications for the LHC physics. In this chapter recent developments within the PYTHIA and HERWIG frameworks 
are reviewed. 


3.1 Multiparton Interactions in Pythia 8 
Contributing author: R. Corke 


3.1.1 MPI in PYTHIA 8 


The original MPI model, first introduced in previous versions of PYTHIA, featured p, ordering, perturbative QCD 
cross sections dampened in the p,; — 0 limit, and a variable impact parameter formalism [8]. These features remain 
in the MPI framework of PYTHIA 8 [56] and have been extended to give a wide range of possible underlying-event 
processes, including all 2 — 2 QCD processes, prompt photon production and others [9]. This newer model was 
developed after the introduction of transverse-momentum-ordered parton showers, opening the way to have a common 
p1 evolution scale for initial-state radiation (ISR), final-state radiation (FSR) and MPI [10]. This common evolution 
is most important for ISR and MPI, which both directly compete for momentum from the beams. 

With such an interleaving, the probability for the it” interaction or shower branching to take place at p1 = p1: is 
given by the combined evolution equation 


dP dPmpt dPisr dPrsr ) 
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where contributions from MPI, ISR and FSR are unitarised by a Sudakov-like factor. The sums for ISR and FSR run 
over all initiator and final-state partons respectively, including those for MPI subsystems, giving full showers from 
these interactions. 

Focusing on just the MPI contribution, the probability for an interaction is given by 


dPmpi _ 1 do oo ( m l doj , ) (3) 
dpı onp dp, p pL ond dp’, ot 


where da/dp , is given by the convolution of PDF factors with the partonic QCD 2 — 2 cross section. This cross 
section is dominated by t-channel gluon exchange, and diverges roughly as dp? /p4 . To avoid this divergence, the 
idea of colour screening is introduced into the model. The concept of a perturbative cross section is based on the 
assumption of free incoming states, which is not the case when partons are confined in colour-singlet hadrons. One 
therefore expects a colour charge to be screened by the presence of nearby anti-charges; that is, if the typical charge 
separation is d, gluons with a transverse wavelength ~ 1/p, > d are no longer able to resolve charges individually, 
leading to a reduced effective coupling. This is introduced by regularising the interaction cross section according to 


dô ag(p) , (Pig + Pi) 
dpł pi (PŽ o +P)? 


, (4) 


where pio (related to 1/d above) is now a free parameter in the model. This parameter is expected to scale with 
energy, and the ansatz is that it does so in a manner similar to the total cross section, an effective power related to 
the Pomeron intercept. The form of the scaling is given by 


E 
E CM 
pio(Ecm) = pi x ( cw) f (5) 


Ow 


where a reference pio is given at some reference energy, and scaled according to Eir 

Up to this point, all parton-parton interactions have been assumed to be independent, such that the probability to 
have n interactions in an event, Pn, is given by Poissonian statistics. This picture is now changed, first by requiring 
that there is at least one interaction, such that there is a physical event, and second by including an impact parameter, 
b. In general, the amount of MPI activity in an event will be directly related to the time-integrated overlap of the 
incoming hadrons during collision, given by 


O(b) = [ou fare Aand aek aes (6) 


after a suitable scale transformation to compensate for the boosted nature of the incoming hadrons. Different matter 
distributions are available, including a single Gaussian, double Gaussian and an intermediate overlap function. While 
requiring at least one interaction results in Pn being narrower than Poissonian, when the impact parameter dependence 
is added, the overall effect is that Pn is broader than Poissonian. 

When the p, evolution has come to an end, colour reconnection is performed. It has been noted, especially by 
Rick Field [57], that a large amount of colour reconnection is necessary to correctly describe data, such as the mean 
pı as a function of charged multiplicity in minimum-bias events. In PYTHIA 8, this is performed by giving each MPI 
subsystem a probability to reconnect with a harder system 


2 
PLlrec 


= y Plrec = R * Dio, (7) 
(PL?ec + Pi) sas 


where R is a user-tunable parameter and p o is the same parameter as in eq. (3. The idea of colour reconnection can 
be motivated by noting that MPI leads to many colour strings that will overlap in physical space. Moving from the 
limit of No —> co to No = 3, it is perhaps not unreasonable to consider these strings to be connected differently due 
to a coincidence of colour, so as to reduce the total string length and thereby the potential energy. With the above 
probability for reconnection, it is easier to reconnect low p} systems, which can be viewed as them having a larger 
spatial extent such that they are more likely to overlap with other colour strings. Currently, however, given the lack 
of a firm theoretical basis, the need for colour reconnection has only been established within the context of specific 
models. 


3.1.2 Rescattering 


A process with a rescattering occurs when an outgoing state from one scattering is allowed to become the incoming 
state of another. This is illustrated schematically in Fig. where (a) shows two independent 2 — 2 processes while 
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(a) (b) 


Fig. 20. (a) Two 2 — 2 scatterings and (b) a 2 > 2 scattering followed by a rescattering 


(b) shows a rescattering process. An estimate for the size of such rescattering effects is given by Paver and Treleani 
[58]. Although rescattering is expected to be a small effect when compared to independent 2 — 2 processes, it should 
be allowed to take place. It would show up in the collective effects of MPI, manifesting itself as changes to multiplicity, 
pı and other distributions, although after a retuning of po and other model parameters, it is likely that its impact is 
significantly reduced. A full implementation of rescattering has been made, allowing the generation of fully hadronic 
final states [59|, and an outline is given below. 

The starting point for the implementation of rescattering is the typical case of small-angle t-channel gluon scatter- 
ing. In this case, the combination of a scattered parton and one of the hadron remnants will closely match one of the 
original incoming hadrons, and the PDF can then be written as 


f(x, Q’) > frescalea(®, Q’) T 5 O(a = Ba) . (8) 


Here, each time a scattering occurs, one parton becomes fixed at a specific £n value, while the remainder is still a 
continuous probability distribution, although rescaled to take into account the momentum taken from it. In this way, 
the original disjoint 2 — 2 MPI are supplemented by single rescatterings, where one parton is taken from the rescaled 
PDF and the other is a delta function, and double rescatterings, where both partons are delta functions. 

Of course, in general, it is not possible to uniquely identify a scattered parton with one of the hadron remnants, 
so an approximate prescription must be used. In particular, we have studied a rapidity based scheme, where at one 
extreme, partons with rapidity y > 0 belong to beam A and those with y < 0 to beam B (“step”) and at the other, 
all partons belong to both beams simultaneously (“sim”). A natural suppression in the amount of single rescattering 
for the simultaneous case means that results do not differ greatly compared to the step prescription. The amount of 
rescattering is also dependent on the amount of underlying activity per event, as the more branchings and scatterings 
in an event, the more partons that are available to rescatter. To study this effect, two different tunes have been studied. 
These are labelled “old” and “new”, where the primary difference is a reduced amount of MPI activity in the new 
tune. 

In Fig. the pı distribution of single rescatterings is shown compared to normal MPI (pp, v/s = 14 TeV) for 
the two tunes. In the ratio plot, double rescattering with the simultaneous prescription and the old tune is also shown 
(with these settings, the effect is largest). As the p, evolution progresses downwards, more and more partons become 
available to rescatter and the rate grows. As expected, however, the rate of single rescattering is small compared to 
normal MPI. As an indicator of the effect of energy on the growth of rescattering, Tab. [1] shows the average number of 
scatterings and rescatterings for different types of event at Tevatron and LHC energies (step option only, old and new 
tunes). Double rescattering is always a very small effect and was neglected in the subsequent hadron-level studies. 

With the full framework implemented, a range of different observables were studied to look for definitive signatures 
of rescattering, including exclusive jet rates, effects on the amount of colour reconnection required to match data, 
enhancement in the pı spectra of final-state hadrons and AR/A¢ distributions. Unfortunately we were unable to find 
any “smoking-gun” signatures. One possible approach to further study where rescattering plays a role would be to 
tune the generator, both with and without rescattering, and to examine differences in the overall fit. 


3.1.3 Tuning prospects 


In this section, tunes of the generator are briefly presented, first for Tevatron data only, and then for an early set of 
LHC data [60]. Parameters relating to final-state showers and hadronisation have been left unchanged, having been 
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Fig. 21. pı distribution of normal MPI and single rescatterings (pp, vs = 14 TeV, old and new tunes). Double rescattering 
with the simultaneous prescription and the old tune is also shown in the ratio plot 


Tevatron LHC 
Min Bias | QCD Jets | Min Bias | QCD Jets 
Scatterings 2.81 5.09 5.19 12.19 
Old Single rescatterings 0.41 1.32 1.03 4.10 
Double rescatterings 0.01 0.04 0.03 0.15 
Scatterings 2.50 3.79 3.40 5.68 
New Single rescatterings 0.24 0.60 0.25 0.66 
Double rescatterings 0.00 0.01 0.00 0.01 


Table 1. Average number of scatterings, single rescatterings and double rescatterings in minimum bias and QCD jet events at 
Tevatron (pp, vs = 1.96 TeV, QCD jet imin = 20 GeV) and LHC (pp, vs = 14.0 TeV, QCD jet fimin = 50 GeV) energies for 
both the old and new tunes 


previously tuned to LEP data. It should be stressed that the parameter space used to make these tunes is somewhat 
limited, and the tunes themselves have been made “by hand” following the principles outlined in [60]. 


Tevatron data 


Two tunes to Tevatron data have been produced, Tune 2C using the CTEQ6L1 PDF set and Tune 2M using MRST 
LO** MOGI]. The MRST LO** PDF set has a relaxed momentum sum rule such that it contains more momentum 
that CTEQ6L1, leading to lower a, and higher p;o parameters. For both tunes, the matter distribution uses an overlap 
function, intermediate between the single Gaussian and the default double Gaussian settings. The balance between 
ISR and MPI activity is primarily tuned by comparisons to pı (Z°) and jet-jet azimuthal angle distributions, both of 
which are driven strongly by ISR. 

In Fig. results from Tune 2C are shown using the Rivet analysis framework for (a) the charged multiplicity at 
ys = 1.8 TeV and (b) the transverse region charged particle density in Rick Field’s leading jet analysis [62157163]. 
Comparisons are made both to data, and to the Pro-Q20 and Perugia 0 tunes of PYTHIA 6 [35/32]. In both cases, 
Tune 2C is in good agreement. 


LHC data 


Tune 4C is a modification of Tune 2C, where MPI and colour reconnection parameters have been changed to give 
a good match to LHC data (with ISR settings left unchanged). Additionally, the diffractive cross section has been 
slightly damped, to better match an early ATLAS study [64]. One notable aspect of this tune is a reduced amount 
of colour reconnection required to match (p,)(No) data in minimum-bias events. This reduction also affects the 
balance between the charged particle density and the scalar py sum density in the transverse region of the underlying 
event, leading to a slightly worse description of this data. 
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Fig. 22. Results from Tune 2C for (a) the charged multiplicity at ys = 1.8 TeV and (b) the transverse region charged particle 
density in Rick Field’s leading jet analysis 
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Fig. 23. Results from Tune 4C compared to ATLAS data at ys = 0.9 and 7 TeV showing (a) the charged multiplicity and (b) 
the transverse region charged particle density 


In Fig. results from Tune 4C are shown against ATLAS data at ys = 0.9 and 7 TeV showing (a) the charged 
multiplicity and (b) the transverse region charged particle density [6516667]. The broad features of the data are 
reproduced, and it is hoped that a more complete tuning will increase agreement further. It is noted that Tunes 2C 
and 2M give too little activity when compared to LHC data and Tune 4C gives too much activity when compared to 
Tevatron data. 


3.1.4 Conclusions 


The MPI model in PYTHIA 8 is an evolution of the original framework introduced in earlier versions of PYTHIA and 
has been well proven in comparisons to experimental data. Of course, it is not the final word in the modelling of 
soft MPI and undoubtedly there are further physics aspects that can be included. One possible extension we have 
implemented is rescattering, outlined above. Although we were unable to find a distinct signature for these processes, 
there is still future scope for the study of effects on an overall generator tune. 

Long awaited LHC data are now being published and will hopefully help to constrain different models and model 
parameters. One current concern is the need for different tunes to describe Tevatron and LHC data, and it remains to 
be seen if there is some region in parameter space where this is possible. So far it is not possible to rule out differences 
due to experimental effects or deficiencies in the model, e.g. related to the energy dependence of pio. We look forward 
to future data that may help to resolve these issues. 


3.2 Multiple Partonic Interactions in HERWIG++ 
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Contributing authors: S. Gieseke, Ch. Rohr, A. Siddmok 


The modelling of underlying events in HERWIG+- is based on the fact that at high enough energies the hard inclusive 
cross section will eventually exceed the total cross section [52]. Generally speaking, we write the hard inclusive cross 
section for partonic 2 > 2 scatters as 


inc min dôi,j 
o(p) =) f dpBfim (e11) @ FH? @ fiyn (ast?) (9) 
ij YPE” t 


which is the usual collinear factorization ansatz. With recent parton distribution functions and lower limits for the 
transverse momentum, p™®i®, in the perturbative regime of a few GeV, Eq. [9] results in values for o”: (s; pe) which 
exceed the Donnachie—Landshoff (DL) parametrization of the total cross section o¢o4. The simplest way out is the 
observation that the proton is a spatially extended object, allowing for independent multiple hard interactions, which 
are strictly all taken into account in the calculation of the inclusive cross section. Therefore, we calculate the average 


number of hard interactions from an eikonal ansatz as 
n(b, s) = A(b; p? )Jo™ (s; pr") . (10) 


The overlap function A(b; u?) describes the spatial overlap of the two colliding hadrons (protons) as a function of the 
impact parameter b. The parameter u? characterizes the size of the proton and is proportional to the squared inverse 
radius. For simplicity, we assume a spatial distribution following the functional form deduced from the electromagnetic 
elastic form factor. But we do allow for a different width parameter u of the distribution, as the colour might be 
distributed differently from the electric charges. The basic ideas for this multiple interaction model follow the model 
in Ref. [II], which in turn introduces a model similar to the one discussed in Ref. [8]. 

The extension to soft scatterings, similar to the model of Ref. [69], is kept as simple as possible. First, the transverse 
momentum of scattered particles is extended to transverse momenta below p@™™. The additional soft contribution to 
the inclusive cross section is also eikonalized, such that we can as well calculate an average number of soft scatters 
from the resulting 0/6, and an overlap function Agog (b) for the soft scattering centers. The functional form Agor(b) is 
assumed to be the same as for the hard scatters, but we allow for a different inverse radius, u2 g. The phase space of 
the soft scatters is determined with a simple Gaussian ansatz for the transverse momentum, of which the parameter 
is determined by the value of o?e. 

The consistency of this model with unitarity is given by fixing the two additional parameters oi™¢, and u2 g from two 
additional constraints. First, we can calculate the total cross section from the eikonal model and fix it to be consistent 
with the DL parametrization. In addition, using the optical theorem, we can calculate the t-slope parameter from the 
eikonal model as well and fix it to a reasonable parametrization. 

After in a first step only the model for hard multiple partonic interactions has been introduced [12], we also studied 
its implications from Tevatron data and total cross section data in a simplified version [70]. Finally, the extension of 
the model to include soft scatters has been implemented in HERWIG++ and is the default underlying event model 
since version 2.3. The consistency of the model predictions with current Tevatron data has been studied in detail in 


Ref. [7]. 


3.2.1 HERWIG+-+ against first LHC data 


Equipped with the good description of the Tevatron data, we can now take a first look at the ATLAS measurements 
made at the 900 GeV and 7 TeV runs at the LHC [65][36]. We anticipate the possibility that the assumptions made in 
order to extend the model into the soft region may be far too simple. Nevertheless, we have been able to accommodate 
the detailed underlying event analyses carried out at the Tevatron. There we have come up with regions in the two- 
dimensional parameter plane of p? and u?, where we obtain a similarly good overall x? for the underlying event data 
and still are consistent with our constraints from the total cross section and the elastic slope parameter. This region 
roughly follows a line. We now had a first look at Minimum Bias data, particularly the relatively simple distribution 
of charged particles in pseudorapidity. 

As a first step, we have varied our model parameters and compared the results against the 900 GeV data. We find 
that the shape of the pseudorapidity distribution in HERWIG+-+ is by far too much peaked in the forward directions. 
In addition, there is not enough freedom in our parameter space to describe (p1 )(Nen). 

A first hint towards the possible improvement of our description of data was found when we varied the probability 
that any of the additional soft scatters gets disconnected in colour space from the rest of the event and the beam 
remnants in particular. The value cD = 1 was used as a default, saying that the soft scatters have always been 
disconnected. Physically this means that there are no colour strings built up between the beam remnants and the 
soft particles produced in the soft underlying event. When they are build up more and more, as we see when we 
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Fig. 24. Pseudorapidity distribution of charged particles in Minimum Bias events at 900 GeV compared to ATLAS data. The 
most sensitive model parameter was the colour disruption probability for soft events. 


vary the parameter towards the other extreme value 0 (always connected), we find that we produce many additional 
soft particles, resulting in an evenly filled plateau in rapidity. Having checked also other parameters, such as parton 
distribution functions and their behaviour at small x values, we found that the effect of the colour disruption parameter 
was most important. Fig. [24] shows the sensitivity to the colour disruption parameter. 

A second hint is given by the unability to describe (p1)(Ncn) which is considered to be very sensitive to non- 
perturbative colour reconnections. So, as final additional modification we have considered a newly implemented model 
for soft colour reconnections in HERWIG+-+. We find that only with the two latter modifications we can give a sensible 
description of minimum bias events. 

A colour reconnection model is a very significant modification of the final state as the production of charged 
particles is affected in its multiplicity as well as in its phase space distribution, once a multiple partonic interaction 
model is used. Hence, before moving on to LHC data, we have checked the new model against data that we previously 
described quite well. First we considered LEP final states and found no sensitivity whatsoever. This was expected as 
the colour structure of the event is well-defined by the perturbative parton shower evolution. The Tevatron underlying 
event analysis, that we were always able to reproduce, was slightly improved, as expected. In particular, the tension 
between pi" and Nen against the leading-jet transverse momentum has been reduced. 

Moving on to the LHC, we have used the ATLAS data with a cut Nen > 6 in order to remove diffractive events 
as we currently have no diffractive model in HERWIG+-+. We have included the variation of the colour disruption 
parameter and the colour reconnection model and tuned our model to the new data. We used data from the run at 
/s = 900 GeV as well as from ys = 7 TeV. The tuned results are shown in Figs. [25] and [26] 

We find that the overall agreement with the data is very good. The expected need for the modification of the colour 
structure in multiple partonic interaction events is in fact given. 


3.2.2 Conclusion 
We have introduced a simple colour reconnection model in HERWIG++ in order to complete the hadronization of 


events with multiple partonic interactions. We find very good agreement with first data on non—diffractive Minimum 
Bias events, measured by ATLAS. The model is included in the recent release HERWIG++ 2.5 [72]. 


4 Theory of Multi-Parton Scattering 


The theoretical investigation of MPI has a long history [73\[74\[75\|4776] and has experienced a renewed interest in 
more recent times [77|(78)/79\/51/[80\81)/821/83\84|[85\[86], driven by the need to understand the hadronic activity at the 
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Fig. 25. Pseudorapidity and charged particle distribution from the ATLAS Nen > 6 analysis at 900 GeV compared to HER- 
wic++. The data points are read off preliminary, but publicly available, ATLAS figures. 


LHC. The phenomenology of multiple parton interactions relies on a rather simple and intuitive cross section formula, 
where multi-parton distributions are multiplied with the cross sections for each individual hard scatter, cf. Chapter [5] 
It is natural to ask whether such a factorization formula can be derived in QCD and to which extent it needs to be 
modified or extended. Sections and [4.2] present two independent efforts to address these questions, respectively 
taking the production of two electroweak gauge bosons or of two jet pairs as examples processes for DPS. 


Accurate predictions of DPS cross sections also require good modelling of double parton distribution functions 
(dPDFs) used in phenomenological studies. Typically, they are constructed from standard single PDFs neglecting 
possible correlations between the longitudinal momenta and transverse positions of the two partons involved. The 
development of the first set of LO dPDFs in the framework assuming factorization between the longitudinal and 
transverse components is described in Section [4.3] of this chapter. However, as discussed there and in Section [4.1] the 
validity of the transverse-longitudinal factorization is being contested. 


Charged particle multiplicity as function of 7 (7 TeV, Nen > 6) Charged particle density (7 TeV, Nu, > 6) 
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Fig. 26. Pseudorapidity and charged particle distribution from the ATLAS Nen > 6 analysis at 7 TeV compared to HERWIG++. 
The data points are read off preliminary, but publicly available, ATLAS figures. 


26 P. Bartalini et al.: Multi-Parton Interactions at the LHC 


This chapter contains also a discussion, in Section of general features of MPI in a probabilistic framework 
employing a functional approach. The probabilistic picture naturally leads to considering the inclusive and exclusive 
cross sections, which are linked by the sum rules. They can be used to obtain information on MPI properties, such as 
two-particle correlations, since inclusive and exclusive cross sections are measured independently. 

The MPI are also very strongly connected to the small-z phenomena of saturation and diffraction. In general, the 
cut diagrams for MPI differ from multiple parton chain diagrams for saturation and diffraction only by the position of 
the cut. In high energy pp scattering unitarity requires the presence of MPI as well as a large diffractive cross section. 
A model to describe diffractive excitations, its relation with the multi-regge formalism and its MC implementation is 
discussed in Section 


4.1 Multiple parton interactions: some theoretical considerations 
Contributing author: M. Diehl 


4.1.1 The basic cross section formula 


As an example for a process to which multiparton interactions contribute we consider the production of two electroweak 
gauge bosons (W, Z or y*) with transverse momenta much smaller than their masses or virtualities. Since the driving 
force for studying multiple interactions is the necessity to describe details of the final state, we keep the cross section 
differential in the transverse boson momenta. For the production of a single gauge boson there is a powerful theory 
involving transverse-momentum dependent parton densities [8718889190], which can to a large extent be generalized 
to the case of multiple hard scattering. 

Figure |27h shows a graph for two-boson production by double parton scattering. This graph can be evaluated 
using the standard hard-scattering approximations, neglecting small momentum components compared with large 
ones. Simple kinematic considerations show that the transverse momenta of the scattering partons are in general not 
equal in the scattering amplitude and its complex conjugate, as indicated in the figure. It is convenient to Fourier 
transform from the transverse-momentum differences r and F to transverse position variables y and y. The constraint 
r +7r = 0 from momentum conservation then turns into y = y for the Fourier conjugate positions, and the cross 
section read] 


do 
Th dx; dx; d?q; 


2 
1 _ _ 
Sg 2 Tl [ee Pk; 5°) (qi — ki — ki) 
Fig. 27} a1 ,42=9,Aq,6q - i=1 

@1,G2=9,Aq,0q 


X G1 ,a;4; (G2) Ô2,azāz (8) fey Faraz (is ki, y) Fanaa (Zi, ki, y), (11) 


where ĉi aia; denotes the hard-scattering cross section for single-boson production. The statistical factor S is 2 if the 
produced bosons are identical and 1 if they are not. Up to power corrections, the momentum fractions of the colliding 
partons are fixed by the measurable momenta as 


zi = q /P = (k} + sr*)/p*, Ti =q /B = (k; + ir )/p ; (12) 


2 In other parts of this chapter, and in Chapter [5] the symbols D or I’ are used to denote the multiple parton distributions 
instead of F used here. Similarly, transverse coordinates are denoted by b, instead of y. 


a b 


Fig. 27. Graphs for the production of two gauge bosons by double (a) or single (b) hard scattering. The dotted line denotes 
the final-state cut. 
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where we have introduced light-cone coordinates v+ = (v? +v3)/v2 for each four-vector v. The definition of the double- 


parton distributions in closely resembles the one for a transverse-momentum dependent single-quark distribution 
BABSIE]. For instance, the distribution of two quarks in the proton with momentum p is given by 


dz; = = 
Foi; ay (Ui, ki, y) = a Ti 0 ži elizi pt k) |2 far 


x (p| G(—§22) Tas 4( $22) A aly + $z) p) 


(13) 


zt =z] =yt= o` 
One can identify k; as the “average” transverse momentum of each quark and y as the “average” transverse distance 
between the quarks, where the “average” refers to the physical scattering amplitude and its conjugate. The result 
thus has an intuitive interpretation: two bosons with transverse momenta q; are produced in the collision of partons 
with average transverse momenta k; and k;. The two collisions occur at an average transverse distance y, which is 
equal to the average transverse distance between the two partons in each colliding proton. At a more formal level, 
F(zxi, ki, y) has the structure of a Wigner distribution in its transverse momentum and position arguments. 
w over the transverse parton momenta, one obtains a collinear two-parton distribution Fy, aa (£i, Y) = 
f Pkid’k Fa as(xi, ki, y), which can be interpreted as the probability for finding two quarks with momentum 
fractions x, and 2x2 at a relative transverse distance y in the proton. These distributions naturally appear if one 
integrates the cross section over the transverse momenta q, and qə of the bosons, 


do 
is dx; dx; 


1 a 2\ a 2 
= S X O1,a1āı (qi) 02,a2ā2 (q2) 
Fig. 27p a1,a2=q, Aq,ôq 
@1,42=4, Agog 


x fev Pai aig (£i Y) Fy, in (Zi; Y). (14) 


This is the simple cross section formula mentioned in the introduction. It was already derived in and underlies 
most phenomenological studies of multiparton interactions. Note that and its generalization (11) to measured 
transverse momenta have an intuitive physical interpretation, but that they arise from the calculation of lowest-order 
Feynman graphs as in figure 27h, using standard approximations. One does not need to appeal to semi-classical 
arguments to obtain these results. 

For each bilinear operator in the matrix element (13) there are three relevant Dirac matrices IG,, 


y= a7": Pag = 375: TÍ = šios with j =1,2, (15) 


which respectively project on unpolarized, longitudinally polarized and transversely polarized quarks. Note that po- 
larized two-parton distributions exist even in an unpolarized proton, where they describe spin correlations between the 
two partons. For small but comparable x; and x2 one may well have sizeable spin correlations between two quarks 
(which are close in phase space for xı ~ x2), even if there is little correlation between the polarizations of a quark 
and the proton (which are far apart in phase space). The relevance of such correlations in multiple interactions was 
pointed out already in [92] but has to our knowledge not been included in phenomenology. 

If parton spin correlations are sizeable, they can have a strong impact on observables. For the production of two 
gauge bosons one can easily see that the product F'aq,Aq F'ag,aq Of longitudinal spin correlations enters the cross 
section with the same weight as the unpolarized term Fy,q F7,q. One also finds that product F%q,5q Fsq,õg of transverse 
spin correlations give rise to a cos(2y) modulation in the angle y between the decay planes of the two bosons and thus 
affects the distribution of final-state particles. 


4.1.2 Power behavior 


It is easy to determine the power behavior of the aoe pochon formula for double hard scattering. The parton-level 
cross sections ĉ behave like 1/Q?, where Q? ~ q? ~ q2 denotes the size w the large squared invariant masses of the 
gauge bosons. We find 
do 1 
MŁ, dx; dz; dq; Fig. [27p Qt A? ? 


where A denotes the size of the transverse momenta qı ~ qə or the scale of non-perturbative interactions, whichever 
is larger. To obtain (16), we have used that the two-parton distributions scale like F ~ 1/A? and that the typical 


(16) 
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Fig. 28. Example graphs for interference terms in fermion number (a) and in quark flavor (b). The blobs indicating two-parton 
distributions are not shown. Labels q and q indicate whether a line is represented by a quark or a conjugate quark field in the 
operator definition. Momenta are assigned as in figure [27] 


transverse distance y between the partons is of order 1/4. The same power behavior as in is obtained for the 
case where both bosons are produced in a single hard scattering, as shown in figure B7þ. Multiple hard interactions 
are therefore not power suppressed as long as one keeps the cross section differential in the transverse momenta of the 
particles produced in the hard parton collisions. 

The situation changes when one integrates over qı and qə. In the double-scattering mechanism both transverse 
momenta are restricted to size A, but for a single hard scattering one has |q, + q| ~ A whereas the individual 
transverse momenta can be as large as Q. Because of this phase space effect one has 


do A? do 1 


n a ~ — ~ =. 17) 
3 - , 2 = ( 
T=: dx; dz; Fig. 27. Qt IG- dx; dx; Fig. 27p Q? 

In the transverse-momentum integrated cross section, multiple hard scattering is therefore a power correction. This is 
required for the validity of the usual factorization formulae, which contain only the single-scattering contribution. 


4.1.3 Additional terms in the cross section 


In the formulae given so far, we have ignored the color structure of the multiparton distributions. The quark lines 
with momentum fraction x; in figure [27h can couple to a color singlet (as in single-parton distributions) but they can 
also couple to a color octet, provided that the lines with momentum fraction x2 are coupled to a color-octet as well. 
In the latter case, the color structure of the operators in is (g Ta, A?q) (q Ta, r°q). Such color-octet distributions 
contribute to the multiple-scattering cross section, as was already pointed out in [76]. They do not have a probability 
interpretation, and little is known about them. The color structure of two-gluon distributions is even more involved. 

There are more multiparton distributions that have no probability interpretation but rather the structure of in- 
terference terms. As shown in figure 28h, the parton with momentum fraction xı can be a quark in the scattering 
amplitude but an antiquark in the conjugate amplitude, provided that the opposite holds for the parton with momen- 
tum fraction x2. Furthermore, one can have interference in the quark flavor, as shown in figure [28p. 

While it is quite straightforward to include these extra contributions in the cross section formula, the additional 
number of two-parton distributions needed to obtain quantitative results is daunting. It is therefore important to have 
some guidance about the size and behavior of these functions. 


4.1.4 High transverse momentum and evolution 


The predictive power of the theory is increased in the kinematic region where the transverse momenta q; are small 
compared with Q but large compared with a typical non-perturbative scale. At least some of the transverse parton 
momenta must then be large as well, and one can evaluate the corresponding parton distributions in terms of a hard 
scattering subprocess at scale |q;| and parton distributions that depend on fewer variables. 

An important example for this are ladder graphs such as the one in figure 29h. Since the partons with momentum 
fraction x; are not connected to those with momentum fraction 72, the momentum mismatch r is small. In position 
space this corresponds to an inter-parton distance y of hadronic size. An important feature of these ladder graphs 
is that their color factors disfavor the color octet distributions mentioned in the previous section. How strong this 
suppression is quantitatively remains to be studied. 
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Fig. 29. (a) Ladder graph for a two-parton distribution, which generates perturbatively large kı and k2 at small r. (b) Graph 
for a quark-antiquark distribution, where the partons with momentum fractions xı and x2 originate from the splitting of a 
single gluon. 


Fig. 30. Graph for the cross section where both two-parton distributions Fa}, ap and F%,,a, (indicated by boxes) involve the 
splitting of one into two partons. 


In the graph of figure [29p, the two partons with momentum fractions x; and x2 originate from the splitting of a 
single parton. Such graphs are relevant for the region of large r, i.e. of short distance y, and also for large kı — kg. 
One finds that they generate strong spin correlations between the partons; the graph in the figure for instance forces 
the helicities of the quark and antiquark to add up to zero. 

Both graphs in figure[29]lead to divergent integrals if one integrates over kı and kz at fixed r. Proper regularization 
of these divergences gives the evolution equation for collinear two-parton distributions; in particular splitting graphs 
as in figure [29p give the famous inhomogeneous term in this equation [93)(94]. Remarkably, one finds however that no 
such inhomogeneous term appears in the evolution of the distributions F(x;, y) at any finite value of y, provided one 
uses the definition with minimal subtraction of divergences. 

The splitting mechanism of figure [29p is responsible for a behavior F'(x;,y) ~ 1/y? at short distances y, which 
renders the integral in the cross section infinite. The same mechanism already gives divergent integrals in the 
differential cross section (11). Furthermore, there is a double counting problem associated with this splitting contri- 
bution: the graph in Anr describes double hard scattering with a 1 — 2 parton splitting in each proton, but 
it also represents a single hard-scattering mechanism, namely two-gluon fusion gg + VV into two bosons via a box 
graph. Therefore, either the cross section formulae and (14), or the definition of two-parton distributions, or 
both must be modified in a way that removes singularities at small y, and a corresponding prescription for calculating 
gg — VV must be given to avoid double counting. This remains a task for future work. More details on the topics 
and results discussed here can be found in [95]. 


4.2 The four jet production at LHC and Tevatron in QCD 
Contributing authors: B. Blok, Yu. Dokshitzer, L. Frankfurt and M. Strikman 


In spite of an extensive theoretical and experimental work, various aspects of the high energy hadronic collisions 
at Tevatron and LHC are still poorly understood. This is especially true for the multijet production which is of a 
paramount importance for the understanding of pQCD dynamics at high energy colliders, and for the search of new 
particles. Here we summarize the first steps of the program to address the topic of MPI starting from the first principles 
of the perturbative QCD which is necessary for an accurate account of the significant disbalance of the momenta of the 
jets (presence of the Sudakov form factors). Among original results presented here are the derivation of the formulae 
in the leading logarithmic approximation for production of 4 jets. Our key finding is that it is possible to isolate the 


30 P. Bartalini et al.: Multi-Parton Interactions at the LHC 


kinematics where the leading twist processes 2 — 4 are not enhanced. This result will allow to improve the reliability 
of the Tevatron studies of the four jet production in the multiparton kinematics and point out directions for the 
corresponding analysis at the LHC. 

Another critical issue is the formulation of the problem in terms of the momentum representation double generalized 
parton distributions and introduction of the mean field approximation for this object. This new formulation is very 
effective for the more detailed studies which are now under way. In addition it brings a link with the original formulation 
in the coordinate space [96//92\[97(98}[99)(1 00)TO1(102)(103)/104)105)(T06)[107], and resolves an issue of the value of the 
strength of the double interaction within this approximation. Previously there was a question whether a conclusion of 
ref. [102\103], that the observed rate is a factor of two larger than the theoretical prediction can be due to uncertainties 
related to many Fourier transforms which were required to convert the HERA data to the experimental number. A 
new formulation, though mathematically equivalent, has completely resolved this issue. This poses serious constraints 
on the Monte Carlo models of pp scattering at collider energies which are not satisfied by many of the current models. 

These issues are of broad interest for both the theorists and experimentalists. 

The standard approach to the multijet production is the QCD improved parton model. It is based on the assump- 
tion that the cross section of a hard hadron—-hadron interaction is calculable in terms of the convolution of parton 
distributions within colliding hadrons with the cross section of a hard two-parton collision. An application of this 
approach to the processes with production of four jets implies that all jets in the event are produced in a hard collision 
of two initial state partons. 

The recent data due to CDF and D0 collaborations do not contradict the dominance of this mechanism in the 
well defined part of the phase space. At the same time these data provide the evidence that there exists a kinematical 
domain where a more complicated mechanism becomes important, namely the double hard interaction of two partons 
in one hadron with two partons in the second hadron. 

Within the parton model picture, the four jets produced this way should pair into two groups such that the 
transverse momenta of two jets in each pair compensate each other. In what follows we refer to this kinematics as 
back-to-back dijet production. We consider the dijets for the case 


is = (ju t+ dae)” < Ja Ihe, Oba < Ja & Ste (18) 


where ô is the total transverse momentum of the dijet and ji the transverse momentum of an individual jet (see 
Fig. (31). The hardness condition 6? >> R~? is implied, with R the characteristic hadron size (non-perturbative scale). 
(The events with disbalances 6? < R~?, give a small contribution both to total and differential cross sections, since 
they are suppressed by Sudakov form factors. (for a detailed detailed discussion of the issue in the review [108] . 
Evidently this nonperturbative contribution is not enhanced in the Leading Logarithmic Approximation). 

Importantly, in this kinematical region the hard scattering of four partons from the wave functions of the colliding 
hadrons remains the dominant source for four-jet production even when the pQCD parton multiplication phenomena 
are taken into account. 

The reason for that is the following. When the two partons from each hadron emerge from the initial state parton 
cascades and then engage into double hard scattering, the resulting differential distribution of the final state jets 
lacks the double back-to-back enhancement factor do œ diz zy which is there in the case of two independent hard 
scatterings. For the two-parton scattering, the characteristic perturbative enhancement do œ 87? results from a 
coherent enhancement of the amplitude due to integration over a large transverse disk, p? ~ 6~? >> jy”. The two 
partons that originate from a perturbative splitting form a relatively compact system in the impact parameter space, 
so that the double hard interaction of such pairs produces only a single perturbative enhancement factor, (613 +624)”, 
which does not favor the back-to-back dijet kinematics (18). The distribution of four jets so produced is much more 
isotropic and can be suppressed by choosing proper kinematical cuts. 

So, the aim of this section is to consider the four-jet production in the hard collisions of four initial state partons. 
We show that the cross section of back-to-back dijet production is calculable in terms of new nonperturbative objects 
— the double-parton Generalized Parton Distributions (DPGPDs). The properties of the DPGPDs can be rigorously 
studied within QCD. In particular, we report here the derivation of the geometric picture for multiple parton collisions 
in the impact parameter space. Up till now, this picture was being used based on a semi-intuitive reasoning 
[98\(99|(TOO/LOTITO2/T03\104\(105)106)(107). 

In the kinematical domain the direct calculation of the light cone Feynman diagrams (momenta of the partons 
in the initial and final states are shown in Fig. using the separation of hard and soft scales shows that the four > 
four cross section for the collisions of hadrons ”a” and ”b” has the form: 


o= f OA fotos fase f dea f azs 


x Da(®1, £2, p2, p2, A) Do (£3, £4, P3, P3, -Ż) 
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Fig. 31. Kinematics of double hard collision - momenta of the colliding partons in in and out states 


d 13 d 24 Op 
7 V dide. (19) 
dt, dt 


Here Da(@1, £2, p2, p2, A) are the new DPGPDs for hadrons ”a” and ”b” defined below. (In the following we will 
consider the case of pp collisions and omit the subscript a. Summing over collisions of various types of partons is 
implied. In practice however we will keep hard scattering of gluons only since it gives the dominant contribution.). 
Remember that the light cone fractions x; are actually fixed by final jet parameters and energy momentum constraints. 

With account of the radiative pQCD effects, in full analogy with the ”DDT formula” for two-body collisions, the 
differential distribution (19) acquires Sudakov form factors [108109] depending on the logarithms of the large ratios 
of scales, j? /8?, and the GPDs become scale dependent: p? ~ 673, på ~ 634. It should be mentioned that the structure 
of the final formula depends on what one actually measures in the experiment — whether energetic single particles 
with large transverse momenta in the final state or ” jets” — and on how the jets are precisely defined. A more detailed 
account of the pQCD effects will be given in a future publication. 


For brevity we will not write explicitly the virtuality scales of the DPGPDs and will use the form: D(x1, x2, Z). 


Note that these distributions depend on the new transverse vector Å that is equal to the difference of the momenta 
of partons from the wave function of the colliding hadron in the amplitude and the amplitude conjugated. Such 
dependence arises because the difference of parton transverse momenta within the parton pair is not conserved. The 
integration limits in x;, f are subject to standard limits determined by kinematic cuts. 

Within the parton model approximation the cross section has the form: 


04 = 0102/7 Rt (20) 


where o; and øz are the cross sections of two independent hard binary parton interactions. The factor mR? charac- 
terizes the transverse area occupied by the partons participating in the hard collision. (In the experimental and 
some of the theoretical papers this factor was denoted as an effective cross section. Our Eq. [21] below shows that such 
wording is not satisfactory since 7R?, does not have the meaning of the interaction cross section.) The data BIA 
indicates that mR? is practically constant in the kinematical range studied at the Tevatron. 


Eq. [19] leads to the general model independent expression for 


(21) 


1 aj PÅ D(e1,22,-A)D(a1, £2, A) 


TR? (2r)? D(a1)D(x2)D(a3) D(a.) 


int 


in terms of DPGPDs. Here D(x;) are the corresponding structure functions. 

The DPGPDs are expressed through the light cone wave functions of the colliding hadrons as follows. Suppose 
that in a four — four process the two partons in the nucleon in the initial state wave function have the transverse 
momenta, Ki, k2. Then in the conjugated wave function they will have the momenta ki + Å, kə — Z. This is because 
only sum of parton transverse momenta but not the difference is conserved. 

The relevant DPGPDs are: 


X f Êk dk 
3 2 — 1 2 2 2 
D (z1, £2, p?, p2, A) z D) (27)? ar~ ky) 
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x (%n (11, k1, £2, ko, ., ki, £i.) 


> > 
x wrt (a1, ki + Aas, = A 5, Berge) + h.c.) 


tan {=n 


x (2r) zi — 1)8(¢ ki). (22) 


i=l i=l 


Note that this distribution is diagonal in the space of all partons except the two partons involved in the collision. Here 
w is the parton wave function normalized to one in a usual way. An appropriate summation over color and Lorentz 
indices is implied. In the case of kinematics 1 >> x1 > x2 we expect only distributions without the spin flip to be 
important. 

Let us stress that it follows from the above formulae that in the impact parameter space these GPDs have a 
probabilistic interpretation. In particular these DPGPD are positively definite in the impact parameter space, cf. 
Eq. Note that in the same way one can introduce the N-particle GPD, Gy, which can be probed in the production 
of N pairs of jets. In this case the first N arguments k; in Eq. [22}are shifted by Z subject to the constraint )7, Zz, = 0. 
So the cross section is proportional to 


i=N 
O2N ~JI T i 
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These GPDs can be easily rewritten in the form of the matrix elements of the operator product. For example: 


D(A) =< N| J d xid zod z3 


x G$, (a1)G%, (w2)G%, (£3)G? , (wa) 
x exp(ipf (a1 — z3) + ips (x2 — z4) 
+ iy (x4 — £3) )|N >, (24) 


calculated at the virtualities p?, p2 at fixed Z. Here we gave an example for the most relevant case of gluons without a 
flip in color and spin spaces. In general a number of distributions can be written, depending on different contractions 
of transverse Lorentz indices and color indices. The classification of the relevant distributions is the same as the 
classification of the quasipartonic operators in ref. [IIO]. Note that the presence of the transverse external parameter 
A does not change the classification, since the corresponding new structures will be strongly suppressed at high 
energies. we wrote the operator expression in the light cone gauge. In arbitrary gauge we shall need Wilson loop W(C) 
connecting points with contracted color indices 
In the approximation of uncorrelated partons it follows from Eq. [22] that 


D(2x1, £2, 7,3, A) = G(21,p}, A)G(z2, p3, A), (25) 


where G(x, Z) are conventional one-particle GPDs. These GPDs can be approximated as Gy (x, Q?, A) = Gn (x, Q?) Fog(A), 
where F2,(A) is the two-gluon form factor of the nucleon extracted from hard exclusive vector meson production (we 
suppress here the dependence of F2, on x) and Gyn (zx, Q?) conventional parton distribution of a nucleon. (Here 

Q? is the virtuality due to the radiation, cf. discussion after Eq. [19]) Thus : 


1 dA m? 
TRZ / (27)? 294) = Sgr i 


Here at the last step we used the dipole fit F2g(4) = 1/(4?/m3} +1)? to the two-gluon form factor. Using the transverse 
gluon radius of the nucleon we obtain 


Ros 2h. re /4 = dFy,(t)/dty=0. (27) 
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This result coincides with the one for the area 7R?,, obtained earlier in using the geometric picture in 
the impact parameter space. That derivation involved taking Fourier transform of the two-gluon form factor and 
calculating a rather complicated six-dimensional integral which could potentially lead to large numerical uncertainties. 
The form of Eq. [27| clearly indicates that numerical uncertainties are small. 


It was emphasized in [102\[103] that the experiments on four-jet production report a smaller value of 7R?,, as 
compared to the one obtained above in the independent particle approximation (though the issue of how well the 
contribution of the 2 — 4 processes was subtracted still remains, cf. discussion in the beginning of this section). It is 
at least a factor of two smaller — that is a four-jet cross section is a factor of two larger — than Eq. [27] gives. (The 
GPDs for sea quarks appear to decrease with A somewhat faster, resulting in a smaller 1/7R?,,, see discussion in 


ant? 
112|.) 
It follows from Eq.[21]that the value of R2, is determined by the range of integration over A. Hence the characteristic 


int 

A in the integral measures the effective distance between the parton pairs (which in principle may differ for different 
flavor combinations). According to the above evaluation within the independent parton approximation the integral for 
1/Rj,, is dominated by small A? ~ 0.1 m2. The contribution of large A is suppressed by the two-gluon form factor of 
a nucleon. This reasoning indicates the important role of inter-parton correlations. In other words, the integral over A 
is effectively cut off by a scale of the nonperturbative correlations. Such correlations naturally arise in nonperturbative 
QCD regime in a number of nucleon models, such as constituent quark model (gluon cloud around constituent quark) 
[L02103], or string model (gluon structure of string). The detailed analysis of the additional correlations due to the 
hard- soft interplay will be reported elsewhere. 


Let us now show that results presented here lead to the intuitive geometric picture in the impact parameter space 
mentioned above [9692//97)/98)[99/TOO/T0N/102\(103)104)[105)T06)/107] . 


The first step is to make transformation into coordinate space i.e., to make Fourier transform from variables k; 
in Eq. [22] to coordinates b;. Performing integration over k; we obtain that transverse coordinates of partons in the 
amplitude and the amplitude conjugated are equal p; = pr. In the calculation we use the fact that upper limit of 
integration over k? is very large compared with the inverse hadron size. Next step is to perform integration over A 
which produces 6(p1 — p2 — p3 + p4) = f d?Bd(pi — ps — B)6(p2 — ps — B). 

The delta functions express the fact that within the accuracy 1/p; where p; is the hard scale, the interactions of 
partons from different nucleons occur at the same point. B is the relative impact parameter of two nucleons. 


The expression for the cross section in the impact parameter space has the form which corresponds to geometry 


of Fig [32] 
o4 = | PBP pepp pDl, £2, P1, P2) 
x D(z3, £4, p3, p4)ô(bı + (B — b3))ô(b2 + (B — b4)) = 
= | @ Ba? pt? pD, 22 pr, 02) 
x D(z3,x4,—B + pı, —B + p2). (28) 
Here the DPGPD in the impact parameter space representation is given by 


D (£1, £2, P1, P2) = 


= 5 Ji [dxd pi] Yn(21, P1, T2, P2, £i, Pis ) 
n=3 


i>3 
x wr (21, P1, £2, P2,---, Li, Pi, OD Tipi). (29) 
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where the delta function expresses the center of mass constraint ea zipi = 0. This is analogous to the case of single 
parton GPDs, see . The functions (#1, 91, £2, P2,.--) are just the Fourier transforms in the impact parameter 
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space of the light cone wave functions and are given by 


=n 


TW @k 
Wn (£1, P1, £2, P2,+--) = Ji Trj xP DL kip) 
i=1 i=1 


x PnalT1, ki, £2, ka, .-)(27) 8 ki). (30) 


Fig. 32. Geometry of two hard collisions in impact parameter picture. 


Thus the GPD defined in Eq. [22]is equivalent to the representation for cross section that indeed corresponds to the 
simple geometrical picture, but instead of a triple integral we now have an integral over one momentum A. Moreover, 
to determinate the cross section we need to know the D(A). The GPD defined in Eq. [22] is useful for calculation of 
many different processes. At the same time the knowledge of the full double GPD is necessary for complete description 
of events with a double jet trigger since the pedestal strongly depends on the impact parameter B [102\(103}. 

Let us stress that this picture is a natural generalization of the correspondence between momentum representation 
and geometric picture for a conventional case of two — two collisions. Indeed in this case it is easy to see that the 
cross section in the momentum representation 


do” a 


a= I Fen) San p) di (31) 


has a simple geometric representation 


do” .. 
a i pid? Bf (es, 1,9?) (22, B - p1, p), (32) 


where f(x, p,p?) = Y+ (a, p, p?)W(a, p, p?) and ~(p, p?) is the Fourier transform of the light cone wave function defined 
above. 

Let us now summarize our results. We have argued that there exists the kinematical domain where the four > 
four hard parton collisions form the dominant mechanism of four-jet production. In this region we calculated the cross 
section, see Eqs. and found that it can be expressed through new two particle GPDs (see Eq. (22h, expressed 
through light cone wave functions. These GPDs depend on a transverse vector A that measures the transverse distance 
within the parton pairs. (Equivalent expressions for these GPDs can be easily given in terms of the operator products.) 
In the impact parameter space we derived the widely used intuitive geometric picture. We argued that the enhancement 
of a four-jet cross section is due to short range correlations in the hadron, as determined by the range of integral of A. 
The contribution of perturbative correlations in the appropriate kinematic domain is suppressed. The detailed study 
of the interplay of the contribution of hard/soft correlations will be reported elsewhere. 

It was argued recently [82\{115], cf. that the cross section can be expressed in terms of two parton distribution 
functions. Our analysis indicates that a more detailed treatment of the QCD evolution effects is necessary. We found 
that it is necessary to introduce the new 2-particle DPGPDs which depend on additional parameter A. The parameter 
A expresses the fact that the difference in transverse components of the parton momenta is not conserved and therefore 
different in |in) and (out| states in the double hard collisions. 


4.3 Double parton distribution functions 
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Contributing author: J. R. Gaunt, W. J. Stirling 


If we make only the assumption that the hard processes A and B may be factorised, then we may write the cross 
section for DPS in very general terms as follows: 


m m a 
(4,3) = > 5 | Talen anb: Qa Qo)tA lEn 21) (02,5) (33) 
Erkal 


xTk (21, £3, b; Qa, Qp)dadr2da\ dr,d”b 


The structure of the DPS cross section formula is similar to that for single parton scattering (SPS), in the sense 
that it too is expressed in terms of parton level cross sections 6 multiplied by parton distributions J’. However, whilst 
the parton level cross sections used in are the same as those used for SPS, the parton distributions of are 
significantly more complex than the single PDFs. They are the generalised two-parton PDFs (or, using the language of 
M6]and Section [4.2} the two-parton GPDs). The quantity I; (£1, £2,b; QA, QB) can be interpreted as the probability 
to find a pair of partons in the proton which have flavours 7 and j, longitudinal momenta xı and z2, and are separated 
in impact parameter by b, when the partons are probed at scales QA and Qp respectively. m is a symmetry factor 
that equals 1 if A = B and 2 otherwise. 

Two further assumptions are often applied to (33). First, it is assumed that the two-parton GPD may be decomposed 
into a (typically flavour independent) transverse piece, and a longitudinal piece: 


Tilt; 73,050 4,08) = F(b) Di (a, 9; Qa, QB) (34) 


In this case the DPS cross section reduces to: 


m ak 
OCA. B) = Srg 5 J devdeade det Di (#1, #2; Qa, QB) Dp (21, 72; Qa, QB) (35) 
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with oe a constant. Second, it is assumed that Dj (x1, £2;QA4, QB) may be written as a product of single PDFs 
(sPDFs): l l 
D} (21, %2;QA, QB) ~ D (21; Qa) D3 (x2; QB) (36) 


The analysis of y +3 jet events performed by the CDF collaboration [3| indicates that approximately holds for 
sea partons at moderately low x. However, it is clear even from elementary considerations that this assumption must 
be violated on some level - for example, the left hand side of must go to zero at the kinematic bound zı + z2 = 1 
whilst the product of sPDFs is finite along this line. One might ask as to whether theory can make any predictions 
on the size of the deviations from (36). 

In 1982, Shelest, Snigirev and Zinovjev derived a ‘double DGLAP’ equation dictating the LO scaling violations of a 
quantity Dij (£1, £2; Q) which we shall refer to as the double PDF or dPDF [94]. Snigirev identified this quantity with 
the factorised longitudinal piece of the two-parton GPD for the case in which the two scales are equal (Q4 = Qg = Q) 
[77]. The equation reads as follows: 
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The first two terms on the right hand side of are associated with changes in the dPDF due to independent 
branching processes — processes in which there are a pair of partons, one of which has the appropriate x and flavour, 
and the other of which splits, either giving rise to the other parton of the appropriate x and flavour, or removing it. 
The most interesting term is the final inhomogeneous term, which represents the increase in the dPDF due to a single 
parton with momentum fraction zı + x2 splitting into a pair with the appropriate x values and flavours. We call this 
the ‘sPDF feed term’ for obvious reasons. The functions P;-,;,;,(z) that appear in this term are known as the 1 > 2 
splitting functions, and may be obtained trivially at LO from the real splitting parts of the usual splitting functions. 
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An important prediction of the double DGLAP equation is that even if the dPDF may be taken to be a product 
of single PDFs at some particular scale Qo, then at any other scale the dPDFs deviate from factorised forms (with 
more deviation the further away one moves from Qo). 

More recently [82], we demonstrated that certain equalities are preserved by the double DGLAP equation, provided 
that they hold at some initial scale Qo. By comparing these equations to an equation in conditional probability, we 
interpreted these equalities as the number and momentum sum rules for the dPDFs, and argued that they should hold 
at the initial scale Qo. The sum rules are: 


Momentum Sum Rule: >| dx1 21D}? (x1, £2; Q) = (1 — 22) D2 (x2; Q) (38) 
—~ Jo 
Jı 
1-22 E Ni DP (22;Q) when jo # jı or Jy 
Number Sum Rule: f dx, Div? (21, 22;Q) = $ (Nj, — UDP (%2;Q) when j2 = jı (39) 
0 


(Ni, + LDP (z2;Q) when jo = jı 
The symbol jiv = jı — jı (jı Æ g), and Nj,, is the number of ‘valence’ jı quarks in the proton. The first sum rule is 
a simple statement of the fact that if one observes a parton with momentum fraction x2 in the proton, the momentum 
fractions of all other partons must add up to 1 — z2. The second rule states that if one observes a parton with flavour 
j in the proton, the number of partons of flavour 7 elsewhere in the proton must be reduced by one (we use the term 
‘number effects’ to describe this simple phenomenon). The sum rules give us some information about the deviations 
of the dPDFs from factorised forms at any scale — certainly, factorised forms do not obey the relations and (39). 
In this contribution, we principally wish to discuss the development of our publicly available set of LO dPDFs — 
the GS09 dPDFs — which have been constructed incorporating pQCD evolution effects and sum rule constraints [82]. 
This discussion may be found in Section However, we should also like to mention some theoretical problems 
that we have recently uncovered in the aforementioned ‘dPDF framework’ for calculating proton-proton cross sections 
— this is covered in Section [4.3.2 


4.3.1 The GS09 dPDFs 


The GS09 dPDF package comprises a grid of dPDF values spanning the ranges 1076 < a, < 1, 107 < x < 1, 
1 GeV? < Q? < 10° GeV’, which is available along with interpolation code from HepForge [117]. It has been obtained 
by constructing inputs that approximately satisfy the sum rules at Qo = 1 GeV, and then numerically evolving these 
inputs up to higher scales according to the d(DGLAP equation. The sPDF set to which we have chosen our dPDF set 
to correspond is (a slightly modified version of) the MSTW2008LO set {118}. 

Given the paucity of experimental data regarding the dPDFs, and in accordance with simple arguments and the 
CDF results, we base our inputs on factorised products of MSTW2008LO sPDFs. However, we modify these basic 
forms in several ways to ensure that the input dPDF's approximately satisfy the sum rules. 

First, all of the dPDFs are multiplied by a factor př (x1, £2) which is designed to take account of phase space 
effects. This factor should ensure the appropriate behaviour of the dPDFs near the kinematic boundary x; + £2 = 1 — 
namely, a smooth decrease to zero. In previous studies, universal phase factors such as (1 — £1 — £2) and (1— 2 — x2)? 
were used — however, with the benefit of knowledge of the momentum sum rules, we can see that neither of these 
options is fully satisfactory. Double PDFs including these factors badly violate the momentum sum rules along the 
lines xı = 0 and z2 = 0. In these regions, a phase factor of approximately 1 would be more sensible (to reflect the 
fact that, ignoring number effects, removing a very low momentum parton does not strongly affect the probability of 
finding any other parton). Indeed, factorised forms satisfy the sum rules perfectly along the lines zı = 0 and x2 = 0. 

It was discovered that the following form for p’! gives rise to inputs which satisfy the momentum sum rules (plus 
appropriate number sum rules) well: 


p” (x1, "2) = (1 — z1 — z2)? (1 — 21)? eO (1 — rg) 2) (40) 


a(t) is 0 if į is a sea parton, and 0.5 if it is a valence parton. The latter two factors on the right hand side are included 
to compensate the decrease of the (1 — 21 — x2)? factor along the lines zı = 0 and x2 = 0. 

We recall that there are only a finite number of valence quarks in the proton, as opposed to an infinite number 
of sea quarks and gluons. Number effects are therefore most significant in the context of valence quarks, and on this 
basis we have chosen to only take account of valence number effects in our inputs. This is done by dividing the uyuy, 
part of any dPDF by two, and completely subtracting the d,d, part. The reasoning behind this is that removing one 
up valence quark essentially halves the probability to find another, whilst there is no chance of finding two valence 
down quarks in the proton. 
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(a) A (b) A 
B B 


Fig. 33. (a): A diagram that apparently contributes to the leading order DPS cross section according to the ‘dPDF framework’ 
described in Section [1] (b): If the dPDF framework of Section{1|is valid, then it should be possible to derive the 1 > 2 splitting 
functions by examining the IR singular parts of this parton-level diagram. 


Finally, we have added extra terms to input distributions whose flavour indices contain jj combinations to take 
account of so-called ‘jj correlations’. These are essentially related to sea parton number effects, although they can 
alternatively be thought of as arising during evolution from some lower scale to Qo via g > jj splittings. It is important 
to include these terms in the equal flavour valence-valence (j,j,) inputs since the ‘jj correlation’ term is much larger 
for these than the quasi-factorised piece at low zx. 

With these adjustments, our dPDF inputs satisfy all sum rules to better than 25% accuracy for x < 0.8 (in the 
normal ‘double human’ basis). 


4.3.2 Is the proton-proton DPS cross section describable in terms of dPDFs? 


Let us consider the calculation of the cross section for a particular (equal-scale) DPS process using the framework 
outlined at the beginning of this section — that is, folding parton level cross sections together with the dPDFs of 
according to (35). Due to the presence of the sPDF feed term in the (DGLAP equation, there will be a contribution 
to the leading order cross section corresponding to diagrams like figure [33[a) in this calculation. For figure [33] a) to 
make a contribution to the leading order DPS cross section, it must be the case that there is a f dk, /k, singularity 
associated with every branching, even at the vertices at which the branches on either side of the diagram split into 
two. Only then will the branches on either side of the diagram contribute a leading (ag log(Q?))" factor to the cross 
section. 

If the dPDF framework is valid, then it should be possible to extract the 1 — 2 splitting functions from a distinct 
part of figure [33{b) that is proportional to (ag log(Q?))?/cer 7, just as it is possible to extract (say) the P,,q splitting 
function from a distinct part of the Drell-Yan diagram with one ISR gluon that is proportional to (ag log(Q?)) [9]. 
However, in [190] (see also [[91]) we demonstrated that for any Standard Model loop with the topology of Batty, there 
is no natural part of the cross section expression corresponding to the loop that has the required structure. We also 
showed more generally that there is no natural part of figure 3a) that is proportional to (ag log(Q?))* /ae ry. 

There is a further, more obvious, problem with the dPDF framework that is actually related to the first problem 
above. In this framework, the same transverse separation profile is effectively assigned to every part of the dPDF, 
with the width of the profile being of the same order as the proton radius. This seems sensible for those parts of the 
dPDF that have arisen as a result of independent branchings. On the other hand, pairs of partons arising as a result 
of perturbative parton splittings do not ‘know’ about the size of the proton — how can it be appropriate to assign an 
effective transverse area of approximately the size of the proton to these pairs? 

There is clearly a flaw in the dPDF framework. The root of the problem seems to be the assumption that the 
two-parton GPD may be decomposed into longitudinal and transverse pieces (as one might have guessed from the 
previous paragraph). Parton pairs arising from perturbative splittings have a singular transverse separation profile, 
which differs significantly from the b profile of partons produced by independent branching [86]. 

The dPDF introduced in Section 1 actually seems to be the integral of the two-parton GPD over b. As a result it 
is directly accessed only in DPS processes in which the two particles probing the proton are uncorrelated in transverse 
space (e.g. the two-nucleon contribution to proton-heavy nucleus DPS [{123)/724]). It is nonetheless interesting to ask 
how ‘wrong’ it is numerically to use plus GS09 dPDFs to calculate DPS cross sections. At present, we can be 
reasonably certain that the contribution associated with the sPDF feed parts of both dPDFs being multiplied together 
should not be included. However, numerically one finds that the accumulated contribution of the sPDF feed term to 
any dPDF is about 10% at low x [82] — therefore this contribution only represents roughly 10% x 10% = 1% of the 
cross section. To definitively answer this question of how wrong the dPDF framework is, we must discover the correct 
framework for calculating proton-proton DPS — we are currently working on this. 
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4.4 MPI: General Features and Consistency Requirements 
Contributing authors: G. Calucci and D. Treleani 


4.4.1 Incoherence and MPI 


MPI cross sections are expressed by the incoherent sum of the contributions due to different numbers of interactions{125] 
[126/59]. A given final state may be however generated by various competing processes, characterized by different 
numbers of partonic collisions and the cross section is the result of diagonal and off-diagonal contributions. Terms 
with different numbers of interactions, giving rise to the same final state, populate the final state phase space in 
a different way and, in the kinematical regions where the contributions to the cross section are similar, important 
interference effects should be expected. 


Fig. 34. a) A diagonal and b) a off diagonal contribution to the MPI cross section. c) an interference and higher order term. 


A diagonal contribution, corresponding to a term with n-partons in the initial state of each interacting hadron, 
is given by the incoherent super-positions of n disconnected parton interactions, localized in n different points in 
transverse space, Fig. 1a). On the other hand, the hard component of the interaction, corresponding to the interference 
between a term with n partons and a term with n’ < n partons, is disconnected and localized in no more than n’ points 
in transverse space, Fig. 1b) [81]. Partons are localized in the hadron by the momenta exchanged in the interaction. 
When partons are localized inside non overlapping regions, much smaller as compared to the hadron size, they are 
only connected one with another through soft exchanges and the picture of independent parallel collisions, which 
characterizes a diagonal contribution to the cross section, is a meaningful one. If, on the contrary, several partons are 
localized by the interaction inside overlapping regions, much smaller as compared to the hadron size, they are allowed 
to interact by exchanging momenta of the size of their virtuality. In an interference term more than two initial state 
partons are localized in the same interaction region|81}. Because of the localization in transverse space, the problem 
of interferences is thus linked to the problem to evaluate the scattering amplitude, at higher orders in the coupling 
constant and including higher twists in the hadron structure, Fig. 1c). One may hence argue that interference terms do 
not represent corrections to the n-pairs of partons scattering inclusive cross section. They rather correct the n’-pairs 
of partons, with n’ < n, scattering inclusive cross section. 

Given the two very different scales in the process, the hadron size and the large momenta exchanged, the hard 
component of the interaction may be disconnected, with the different hard parts linked only through soft exchanges 
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and localized in different regions in transverse space. Different MPI terms are hence conveniently understood as the 
contributions to the final state due to the different disconnected regions, where the hard component of the interaction 
is localized. In the simplest case, in each different disconnected region the interaction may be evaluated at the lowest 
order in the coupling constant. In each single partonic collision all transverse momenta balance and a MPI process 
contributes to the cross section by generating different groups of final state partons where the large transverse momenta 
are compensated separately. 

When MPI are understood in the topological sense described above, different MPI terms, corresponding to different 
localizations in transverse space, do not interfere and the final cross section is obtained by the simple superposition 
of the cross sections, due to the contributions of the different topological configurations of the hard component of the 
interaction. MPI hence add incoherently. A remarkable consequence is that MPI allow a probabilistic description. 


4.4.2 The Probabilistic Picture of MPI: A Functional Approach, Cancellation of Unitarity Corrections 


As Multiple Parton Interactions add incoherently, the problem may be discussed within a probabilistic framework|5]] [80]. 
A functional approach is most general. One may start by introducing the probabilities W,,, to find the hadron in con- 
figurations with n-partons with coordinated u1... Un, ui = (bi, £i), where b; are the transverse parton coordinates and 
x; the fractional momenta, and the multi-parton generating functional, Z: 


Wii), 20 =S 4 | He) Hun) Wa (us + tn de dy (41) 


Probability conservation implies the normalization condition Z[1] = 1, while the probabilities of the various configu- 
rations, W,,, are the coefficients of the expansion of Z[J] at J = 0. The coefficients of the expansion of Z[J] at J = 1 
give the many body densities, i.e. the inclusive distributions: 


622 
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Correlations, which describe how much the distribution deviates from a Poissonian, are obtained by the expansion of 
the logarithm of the generating functional, F[J] = InZ[J], at J = 1: 


F|J] = [Peo —1jdu+ SO = f -< Un) [F(u1) 1]... 
_ | [F(in) — dune dtia (43) 


One has F[1] = 0 and, in the case C, = 0,n > 2, the multi-parton distribution is a Poissonian. 
Given the multiparton distributions W,,, one may express the hard cross section in a functional form: 
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n,m 


x {1 — i I [1 — ĉ; z(u, u’)] } ] | duav' 


i=1 j=1 


Zp|J'] 


(44) 


J=J'=0 


here £ is the impact parameter between the two interacting hadrons A and B and ĝ; j is the probability for the parton 
i (of A) to have an hard interaction with the parton j (of B). 

The hard cross section is obtained by summing all contributions due to all different hadronic configurations (the sums 
over n and m). For each pair of values n and m, one has a contribution to Chara when at least one hard interaction 
takes place. Given n and m, the probability to have at least one hard interaction is represented by the term in curly 
brackets. 

The cross section is analogous to the expression of the inelastic nucleus-nucleus cross section in the Glauber model{127] 
and takes into account both disconnected interactions (which imply n = m) and rescatterings (when n 4 m). In 
the present case one is interested in disconnected interactions. One may hence simplify the expression in Eq.(4) by 
neglecting all rescatterings. To that purpose the term in curly brackets is replaced with the following expression: 
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where, in the second line, all repeated indices (which correspond to rescatterings) have been suppressed. In particular 


only the first term of the expansion of the logarithm in the first line in Eq.(5) needs to be taken into account. The 
symmetry of the integrand with respect to the indices allows obtaining a compact expression of the hard cross section: 


Cnara(B) = exp(ð) -exp(ð") f1 — exp(-d-6-8’)] Za1Za17] 
Š [1 - exp(- 0.6. az alJ}Zs[J'] 


J=J'=0 
(46) 


J=J'=1 


where all convolutions are understood. 
Eq. (6) includes all MPI, which are identified with the disconnected collisions. The hard cross section Ghara is easily 
expressed as a sum of MPI: 
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where 6(212',) in the last line of Eq.(8) is the parton-parton cross section integrated with a cutoff. Given the localization 
of the interactions in transverse space, the parton-parton interaction probability has been treated as a d6-function of 
the transverse coordinates, namely ô(u, u’) = G(a, 2')d(b — b’). 

The average (N)onara is hece equal to og, the single scattering inclusive cross section of the QCD parton model. 
In an analogous way one obtains 


N(N-1 1 A 5 
A 6 sara B) T ay | Parez biba) (ora Jor) 
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where øp is the double parton scattering inclusive cross section. The relation is easily extended to any number of 
MPI. One obtains 


a1...2K301...bK)6(2121)...6(2K 2) 
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K 
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which proves that, when rescatterings are neglected, for any choice of multiparton distributions, the inclusive cross 
sections are given by the moments of the distribution in the number of collisions. 
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4.4.3 Inclusive and Exclusive Cross Sections, Sum Rules 


In proton-proton collisions, the inclusive cross sections are thus the moments of the distribution in the number of MPI. 
The most basic information on the distribution in the number of collisions, the average number, is hence given by 
the single scattering inclusive cross section of the QCD parton model. Analogously the Ath scattering inclusive cross 
section gives the Kth moment of the distribution in the number of collisions and is related directly to the K-partons 
distribution of the hadron structure. 

A way alternative to the set of moments, to provide the whole information of the distribution, is represented by 
the set of the different terms of the probability distribution of multiple collisions. Correspondingly, in addition to the 
set of the inclusive cross sections, one may consider the set of the ”exclusive” cross sections, where one selects the 
events where only a given number of collisions are present|5I]. The cross sections called now ”exclusive” are in fact 
partially inclusive cross sections, since one sums over all large p, partons outside a given phase space interval and over 
all soft fragments. 

Interestingly, in its study of double parton collisions, the CDF experiment[3] did not measure the double parton 
scattering inclusive cross section. The events selected where in fact only those containing just double parton collisions, 
while all events with triple scatterings (about 17% of the sample of all events with double parton scatterings) where 
removed. The resulting quantity measured by CDF is hence different with respect to the inclusive cross sections usually 
discussed in large momentum transfer physics. In fact CDF measured the double parton scattering ” exclusive” cross 
section. 

One should emphasize that the ” exclusive” cross sections are not given by the usual QDC-parton model expression 
of large p, processes. While the inclusive cross sections are in fact linked directly to the multi-parton structure of the 
hadron, the link of the ” exclusive” cross sections with the hadron structure is much more complex. The requirement 
of having only events with a given number of hard collisions implies in fact that the corresponding cross section (being 
proportional to the probability of not having any further hard interaction) depends, at least in principle, on the whole 
series of multiple hard collisions and hence on an infinite non-perturbative input. One has 


ae X N(N-1)...(N-K +1). 
Thara = >, ON, OK = 5 ( zm den (51) 
N=K 


where Gy is the ”exclusive” cross section of N partonic collisions, while ox is the inclusive cross section of K partonic 
collisions. Notice that Eq. Bilis a set of sum rules that relate inclusive and ”exclusive” MPI cross sections. Notice also 
that MPI add incoherently in the final cross section, leading to a different probabilistic picture of the process in each 
phase space interval, in such a way that one may associate a different probability distribution of MPI to each different 
phase space choice of observing the final state. In a given phase space window only a small number of MPI might give 
a sizable contribution to Chara. If, as an example, the relevant contributions to Chara are at most from triple collisions, 
one may write 


0, = os —20p + 30r 
Õ2 = 0p— 307 
03 = oy (52) 


which shows how to express the ”exclusive” cross sections in terms of well defined elements of the non perturbative 
hadron structure, in the phase space window under consideration. 

Any final state phase space window identifies an interval in momentum transfer and in fractional momenta, which 
represents the domain of definition of the corresponding probability distribution of MPI. The same interval in momen- 
tum transfer and fractional momenta represents also the integration domain of the integrated terms, which appear in 
the exclusive” differential cross sections. Integrated terms appear in the ”exclusive” differential cross sections because 
of the normalization of the probability distribution. The limits of the integrated terms are thus fixed unambiguously 
by normalization and coincide with the kinematical limits adopted to select the final state. 
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In Eq.(13) the explicit expressions of the single scattering inclusive, og(u, u’), and ” exclusive”, 1 (u, u’), differential 
cross sections are given as a function of the coordinates u, u’ of the observed partons. The expression of the inclusive 
cross section holds at all orders in 6, while the ” exclusive” cross section is at the order 6?. 


Single scattering differential “exclusive” D,, O 
cross sections at order 6? including the 


two-body parton correlations: 
C4 B 0 


average number of 
partons 


correlation in a 


(at all orders in Ĉ) 


A B A pair 
Oo — 
O (ô ) integrated cross 
5! == section 
Oman) differential cross 
A B section 
Exclusive Q 2 
diff. cross 7 i 
section Q © 
A B A B 


Fig. 35. Single scattering differential inclusive and ” exclusive” cross sections including the two-body parton correlations. The 


* exclusive” cross section is at order 67. 


By restricting the phase space interval of the observed final state, a, is therefore well expressed by the term linear 
in ô. In that limit the probability of interaction is well approximated by the average of the distribution, while 2 is 
negligibly small and the single parton scattering ”exclusive” cross section is well represented by the single scattering 
expression of the simple QCD parton model og. When the phase space volume is increased, the single scattering 
” exclusive” cross section becomes increasingly different from the single scattering of the simple QCD parton model 
and the difference between og and 0; allows a direct measure of the importance of correlations, Fig.2. 

Notice that the components of the hadron structure, namely the terms D and C, are well defined quantities, as 
they do not mix when the kinematical limits adopted to select the final state are changed. The effect of modifying 
the kinematical limits adopted to select the final state is in fact only to change the integration domain of each term, 
namely to probe the multi-parton structure of the hadron in different domains in x and Q? [BI]. 

Inclusive and ”exclusive” cross sections result from independent measurements. By checking the sum rules, com- 
paring the measured integrated inclusive cross section with the sum of the measured integrated ” exclusive” cross 
sections, taken with the proper multiplicity factors, up to a given order in the number of collisions, one has hence a 
direct indication of the importance of terms with larger numbers of partonic collisions in a given phase space interval. 
In a phase space interval where, as an example, only single and double collisions give sizable contributions, one may 
thus obtain information on the effect of two body parton correlations by looking at the difference between the single 
scattering inclusive and ”exclusive” differential cross sections: 


dog da, = dos Os (54) 
dydp; dydp; dydp; Ceff 


By comparing the behavior, as a function of fractional momenta, of the difference in the left hand side of Eq.(14) with 
the right hand side one obtains information on the value of the effective cross section and on the dependence of oes 
on y and pr. 
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4.4.4 Concluding Summary 


Given the extended nature of the interacting systems, the hard component of a hadronic interaction is, in general, 
disconnected. The different terms of a MPI process correspond to the different regions where the hard component of the 
interaction is localized. In a single partonic collision all transverse momenta balance and a MPI process contributes 
to the cross section by generating different groups of final state partons where the large transverse momenta are 
compensated separately. A consequence is that terms with different numbers of hard interactions add incoherently and 
the MPI cross section is obtained by the superposition of the cross sections, due to the contributions of the different 
topologies of the hard component of the interaction. 

The physical picture of MPI is hence probabilistic, which leads naturally to consider two different sets of cross 
sections, the inclusive cross sections, given by the moments of the probability distribution of MPI, and the ” exclusive” 
cross sections, given by the different terms of the probability distribution. The inclusive cross sections depend on all the 
different terms of the distribution in multiplicity of partonic interactions, each counted with a different multiplicity 
factor, while each ”exclusive” cross section corresponds to a given term of the distribution. A consequence is that 
inclusive and ”exclusive” cross sections, which result from independent measurements, are linked by sum rules, which 
are implied by the relations expressing the moments of a probability distribution. 

Of course different final state phase space windows lead to different probability distributions of MPI and the number 
of MPI can be controlled by adjusting the final state phase space interval. The sum rules are thus saturated with a 
different number of terms in each different final state phase space interval. The number of terms needed to saturate 
the sum rules in the phase space interval under consideration provides a quantitative measure of the importance of 
terms with larger numbers of partonic interactions. Given the number of terms needed to saturate the sum rules, the 
non-perturbative input to the ”exclusive” cross sections is given explicitly in terms of well defined properties of the 
hadron structure, which hence allow to evaluate unambiguously also the ”exclusive” cross sections in perturbative 


QCD. 


4.5 Multiple interactions, diffraction, and the BFKL pomeron 
Contributing author: G. Gustafson 


In high energy pp scattering the cross section for minijet production becomes very large, and unitarity implies that 
multiple interactions, saturation, and diffraction become important. Effects of saturation and multiple interactions are 
most easily described in impact parameter space, as parton rescattering is represented by a convolution in transverse 
momentum space, which corresponds to a simple multiplication in transverse coordinate space. 

The proton has an internal substructure, which may be excited in a diffractive scattering process, and diffractive 
excitation represents large fractions of the cross section in pp collisions or DIS. In the Good—Walker formalism [128] 
diffractive excitation is described by the fluctuations in the scattering amplitude. In most analyses of pp collisions this 
mechanism is used only for low mass excitation, while high mass excitation is described by a triple-regge formula [129] 
130], where regge trajectories and couplings are fitted to experimental data (for recent analyses see e.g. Refs. [131] 
M3233). 

The proton substructure is represented by a parton cascade, which at high energies is described by BFKL evolution. 
The fluctuations in this evolution are known to be very large [134]. An analysis of these fluctuations, within the Lund 
Dipole Cascade model, is able to reproduce the experimental cross sections for diffractive excitation in pp collisions 
or DIS [135136]. This implies that the effective pomeron couplings in the multi-pomeron formalism can be estimated 
without any new free parameters. 

It should also be noticed that the classification of diffractive events varies between different formalisms, and cannot 
be uniquely defined. Therefore it is recommended to study gap events rather than diffractive events. 


4.5.1 The eikonal approximation and the Good—Walker formalism 


As mentioned in the introduction, diffraction, saturation, and multiple interactions are most easily described in impact 
parameter space. If the scattering is driven by absorption into inelastic states 7, with weights 2f;, the elastic amplitude 
is given by 

T=1-e", with F=} fi (55) 


For a structureless projectile we find: 
do tor /d7b = (2T), 


ca /2b = T), (56) 
Tinel / db = (1 = eT fi) ~N Otot — Oel. 
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If the projectile has an internal structure, the mass eigenstates Y% can differ from the eigenstates of diffraction 
®,, which have eigenvalues Ta. With the notation Wp = >>), ckn@n (with Win = Y1) the elastic amplitude is given 
by (|T|) = S>c4,Tn = (T), while the amplitude for diffractive transition to mass eigenstate YW is given by 
(Y%|T |W) = Jn CkenTncin. The corresponding cross sections become 


doe /@b = (X Rpa}? = (T) (57) 
dowis¢/d°b = X(W |T We) (PT) = (T°). (58) 
k 


The diffractive cross section here includes elastic scattering. Subtracting this gives the cross section for diffractive 
excitation, which is determined by the fluctuations in the scattering process: 


dodit f ex/ db = doaif f > da el = (T?) = Py (59) 


4.5.2 Dipole cascade models 


Mueller’s dipole cascade model [137\(138)[139] is a formulation of BFKL evolution in transverse coordinate space. Gluon 
radiation from the color charge in a parent quark or gluon is screened by the accompanying anticharge in the color 
dipole. This suppresses emissions at large transverse separation, which corresponds to the suppression of small kq in 
BFKL. For a dipole (x,y) the probability per unit rapidity (Y) for emission of a gluon at transverse position z is 
given by 

dP a Bz (x — y)? 305 


= ives 
dY 2r (x-z}(z- y)? WH aT og 


(60) 


The dipole is split into two dipoles, which (in the large Ne limit) emit new gluons independently. The result is a 
cascade, where the number of dipoles grows exponentially with Y. 

When two cascades collide, a pair of dipoles with coordinates (x;,y;) and (xj, yj) can interact via gluon exchange 
with the probability 2fij, where 


ursis (Eten o 


Summing over all dipoles in the cascades then reproduces the LL BFKL result. The elastic scattering amplitude is 
given by T = 1 — exp(— Ì_ fij), and the cross sections are given by Eqs. 

The Lund cascade model is a generalization of Mueller’s model, which includes: 

— NLL BFKL effects 

— Nonlinear effects in the evolution 

— Confinement effects 

For an incoming virtual photon splitting in a q@ pair, the initial state wavefunction is determined by perturbative 
QCD. For an incoming proton we make an ansatz in form of an equilateral triangle of dipoles. After evolution the 
result is rather insensitive to the exact form of the initial state. The model is also implemented in a MC program 
DIPSY. The model reproduces successfully the total and (quasi)elastic cross sections for DIS and pp scattering. 


4.5.3 Fluctuations and diffractive excitation 


The fluctuations in the evolution are large, and the model can also describe diffractive excitation within the Good- 
Walker formalism, without new parameters beyond those adjusted to the total and elastic cross sections [135]. This is 
similar in spirit to the early analysis by Miettinen and Pumplin [143]. In DIS saturation effects are not very important, 
while in pp collisions saturation effects strongly suppress the fluctuations, and thus the cross section for diffractive 
excitation. 

y*p collisions 

The distribution in the non-saturated scattering amplitude, F, is shown in Fig. [36] for different impact parameters. 
The distribution can be approximately described by a power qe x AF? (with a cutoff for small F-values), which is 
illustrated by the straight lines in the figure. The width of this distribution is rather large, and the approximation gives 
the ratio dadis f.ex./dotot % 1—1/ 2?-P In the simulations the power p is rather independent of the impact parameter, 
and therefore this result is also valid for the integrated cross sections. This gives oaip¢/Otot ~ 0.13 at Q? = 50 GeV?, 


decreasing for larger Q?, but fairly insensitive to the energy W. 
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Fig. 36. Distribution in the one-pomeron amplitude F in DIS for Q? = 14 GeV? and W = 220 GeV. The photon is here 
represented by a dipole with size r = 1/Q. The impact parameter b is measured in Gev~?. 
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Fig. 37. Distribution in the one-pomeron amplitude (F', left) and the uniterized amplitude (T, right) in pp collisions at 2 TeV. 
Notation as in Fig. 1. 


pp collisions 

In pp scattering the Born amplitude is large, and therefore unitarity effects are important. Figure [87] shows both 
the Born amplitude and the unitarized amplitude at 2 TeV for different b-values. We see that the width of the Born 
amplitude is large, and without unitarization the fraction of diffractive excitation would be similar to that for y*p for 
lower Q?-values. (The smooth lines are fits of the form AF?e~*".) 

However, the unitarized amplitude is limited by 1, and the width, and therefore the diffractive excitation, is very 
much reduced. This is in particular the case for central collisions where the amplitude approaches the black disc limit, 
as seen in the right panel in Fig. This result corresponds to the effect of enhanced diagrams in the conventional 
triple-regge approach. The impact parameter profile is shown in Fig. We see that the cross section for diffractive 
excitation is largest in a ring with radius b ~ 1fm. This result also implies that factorization is not satisfied when 
comparing diffractive excitation in DIS and pp scattering. 


4.5.4 Comparison with multi-regge analyses 


It is also interesting to compare the results from the Good—Walker analysis with the multi-regge formalism. To this end 
we study the contribution from the bare pomeron, meaning the one-pomeron amplitude without contributions from 
saturation, enhanced diagrams or gap survival form factors. When s, MZ, and s/M@ are all large, pomeron exchange 
should dominate. The results of the MC for the total, elastic, and single diffractive cross sections are shown by the 
crosses in the right panel in Fig. We note that the results are very well reproduced by a triple-regge expression 
with a single pomeron pole, with the parameters 


a(0)=1.21, a’ =0.2GeV~*, gpp(t) = (5.6 GeV™t) et, g3p(t) = 0.31 GeV’, (62) 


which is shown by the straight lines. Also the t-dependence of diffractive excitation is well reproduced. 
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Fig. 38. Left: Impact parameter distributions for (T) = (dotor/d?b)/2, (TY? = doei/d?b, and Vr = doaisex/d?b in pp collisions 
at W = 2 TeV. b is in units of GeV~'. Right: The total, elastic, and single diffractive cross sections in the one-pomeron 
approximation. The crosses are model calculations and the lines are from a tuned triple-regge parametrisation. 


These results can also be compared with multi-regge analyses, where e.g. Ryskin et al. obtain a(0) = 1.3, 
a’ < 0.05GeV~?, Kaidalov et al. [132] find a(0) = 1.12, a’ = 0.22GeV~?, while Gotsman et al. [133] find a(0) = 
1.335, a’ = 0.01 GeV~?. Thus our values are somewhere in between. We also note that the Good—Walker results are 
reproduced by a single pomeron pole, i.e. not by a cut as expected in LL BFKL, or a series of poles as obtained 
with a running coupling [144]. Also the triple-regge coupling g3p is approximately constant, while in LL BFKL it is 


proportional to ~ 1/,/|t| [T38)/745). 


4.5.5 Can diffraction be uniquely defined? 


Multipomeron diagrams are included in the dipole picture, with fixed multi-pomeron couplings. However, all events 
with no gap are classified as inelastic. This can be compared to the formalism in Ref. [131], in which a large cross 
section for overlapping double diffraction is obtained. 

To compare predictions from different models with each other and with experiments, we need a unique definition 
of diffraction. One attempt might be events with two separate color singlet systems, containing the original valence 
quarks. This could be obtained by exchange of two (or more) gluons forming a color singlet. If such a state is obtained 
in a calculation in perturbative QCD, a gap could, however, be filled by final state radiation or nonperturbative 
strings, or else a gap could be formed by color reconnection. Thus diffractive events cannot be uniquely defined by 
perturbative QCD. The definition varies between different schemes, and for a specified event the diffractive capacity 
is not an observable. The solution must be: study observables, meaning events with a gap. 


4.5.6 Summary 


— In high energy pp scattering unitarity implies a high probability for multiple interactions, and a large diffractive 
cross section. In the Good—Walker formalism elastic scattering is given by the average of the scattering amplitude, 
while the fluctuations determine the diffractive excitation. 

— Parton cascades fill the whole rapidity range between projectile and target. The fluctuations in BFKL evolution are 
large, and within the Lund Dipole Cascade model they can describe diffractive excitation within the Good—Walker 
formalism, to both low and high masses. 

— When the interaction approaches the black disc limit in central pp collisions, fluctuations and diffractive excitation 
is suppressed. This leads to factorization breaking in comparisons of DIS and pp scattering. 

— The result of the model calculations corresponds to a bare pomeron, which is a simple pole, and an almost constant 
triple-pomeron coupling. 

— Diffractive excitation is scheme dependent, and cannot be uniquely defined. Study gap events. 


5 Phenomenology of Multi-parton Interactions 


Due to its high collision energy and luminosity, the LHC provides a valuable opportunity to observe multiple parton 
hard-scatterings, in particular many DPS processes. Theoretical investigations of double parton scattering have a long 
history, with a large number of studies evaluating the DPS contribution to high energy processes [96][58)/146)/14 7/148] 
(149/71 50\15 1752\/153)154)1155)[156) . 
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Compared with the DPS processes already observed, namely final states involving 4 jets (at the AFS collaboration 
at the CERN ISR [I]), and y + 3 jets (at the CDF [3] and the DO [4] collaborations at the Fermilab Tevatron), processes 
at the LHC involve different scales and initial state partons, hence providing complementary information on DPS. 
Moreover, large contributions from double parton scattering could, for example, result in a larger rate for multi-jet 
production than otherwise predicted, and produce relevant backgrounds in searches for signals of new phenomena. It is 
thus important to know empirically how large the double parton contribution may be and what dependence on relevant 
kinematic variables it has. Specifically, the differences between final states produced in SPS and in DPS processes need 
to be studied in order to separate between these processes and gain more detailed experimental information on DPS. 
In addition to its role in general LHC phenomenology, the DPS measurements will have an impact on the development 
of partonic models of hadrons, since the effective cross section for double parton scattering measures the size in impact 
parameter space of the incident hadron’s partonic hard core. 

The DPS phenomenology is based on the general expression for the cross section ota. B) cf. Section 4.3] 


ota By) = T 5 J exrdcade dye? bP (En, 225b tista) Fae hsbs tr: ta) 01, 04 )98 (ea, 0) (63) 
ijkl 


The I; (£1, £2, b; t1, t2) represent double parton distributions. They may be loosely interpreted as the inclusive proba- 
bility distributions to find a parton i(j) with longitudinal momentum fraction x1(x2) at scale tı = In(Q?)(t2 = In(Q3)) 
in the proton, with the two partons separated by a transverse distance b. The scale t;(t2) is given by the characteristic 
scale of subprocess A(B). The quantity m is a symmetry factor that equals 1 if A = B and 2 otherwise. Separating 
the transverse part, T; j (£1, 22, b; ti, t2) = Di j(£1, V2, t1, tz) x F(b), Eq. |63| becomes 


DPS m 1 1 l pi aA 1 aB 1 
(A,B) = zg » J exrdrade da Diy (21, 22, 1, ta) Dealt th ta, ta)OA( a1 24th) Ft (22: 2ta), 
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If one makes the further assumptions that double parton distributions reduce to the product of two independent one 
parton distributions, D; ; = D; x D;, the DPS cross section can be expressed in the simple form 


Teff 


DPS _ MOACB 
o =— : (65) 
(5) 2 Gest 
In this section the phenomenological studies for the production of two jets in association with a bb pair (in section 
5-1), same-sign W pair (in section|5.2) and Z+ jets (in section |5.3) are reviewed. 


5.1 Dynamical Characteristics of Double Parton Scattering 
Contributing author: E. Berger 


The concept of a DPS process consisting of two short-distance subprocesses occuring in a given hadronic interaction, 
with two initial partons being active from each of the incident protons, is shown for illustrative purposes in Fig. It 
may be contrasted with conventional single parton scattering (SPS) in Fig/40| in which one short-distance subprocess 
occurs, with one parton active from each initial hadron. Both produce the same 4 parton final state. Our aims 
in Ref. are to address whether double parton scattering can be shown to exist as a discernible contribution 
in well defined and accessible final states, and to establish the characteristics features that allow its measurement. 
We show that double parton scattering produces an enhancement of events in regions of phase space in which the 
“background” from single parton scattering is relatively small. If such enhancements are observed experimentally, with 
the kinematic dependence we predict, then we will have a direct empirical means to measure the size of the double 
parton contribution. 

From the perspective of sensible rates and experimental tagging, a good process to examine should be the 4 parton 
final state in which there are 2 hadronic jets plus a b quark and a b antiquark, viz. b b jı j2. If the final state arises 
from double parton scattering, then it is plausible that one subprocess produces the b b system and another subprocess 
produces the two jets. There are, of course, many single parton scattering (2 to 4 parton) subprocesses that can 
result in the b b jı jz final state, and we identify kinematic distributions that show notable separations of the two 
contributions. 

The state-of-the-art of calculations of single parton scattering is well developed whereas the phenomenology of 
double parton scattering is less advanced. For pp — bbj1j2X, assuming that the two subprocesses A(i j — b b) and 
B(kl— jı j2) in Fig. [B9]are weakly correlated, and that kinematic and dynamic correlations between the two partons 
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Fig. 39. Sketch of a double-parton process in which the active partons are i and k from one proton and j and l from the second 
proton. The two hard scattering subprocess are A(i j > a b) and B(k l > c d). 


Fig. 40. Sketch of a single-parton process in which the active partons are į from one proton and j from the second proton. The 
hard scattering subprocess is A(i j > a b c d). 


from each hadron may be safely neglected, we employ the common heuristic expression for the DPS differential cross 
section 


= d SPS bbX d SPS ae 
doP?S (pp = thf joX) = oe > DANE ep A) (66) 


The numerator is a product of single parton scattering cross sections. In the denominator, there is a term cep with 
the dimensions of a cross section. Given that one hard-scatter has taken place, ce measures the size of the partonic 
core in which the flux of accompanying short-distance partons is confined. Collider data [4] yield values in the range 
Cert ~ 12 mb. We use this value for the estimates we make, but we emphasize that the goal should be determine its 
value at LHC energies. 

In Ref. [157], we present the details of our calculation of the double parton and the single parton contributions 
to pp —> b b jı j2 X. We perform full event simulations at the parton level and apply a series of cuts to emulate 
experimental analyses. We also treat the double parton and the single parton contributions to 4 jet production, again 
finding that good separation is possible despite the combinatorial uncertainty in the pairing of jets. 


5.1.1 Distinguishing variables 


Correlations in the final state are predicted to be quite different between the double parton and the single parton 
contributions. For example, we examine the distribution of events as function of the angle @ between the planes 
defined by the bb system and by the jj system. If the two scattering processes ij — bb and kl > jj which produce the 
DPS final state are truly independent, one would expect to see a flat distribution in the angle ®. By contrast, many 
diagrams, including some with non-trivial spin correlations, contribute to the 2 parton to 4 parton final state in SPS 
ij — bbjetjet, and one would expect some correlation between the two planes. In Fig. we display the number of 
events as a function of the angle between the two planes. There is an evident correlation between the two planes in 
SPS, while the distribution is flat in DPS, indicative that the two planes are uncorrelated. 

Another interesting difference between DPS and SPS is the behavior of event rates as a function of transverse 
momentum. As an example of this, in Fig. we show the transverse momentum distribution for the leading jet 
(either a b or light j) for both DPS and SPS. SPS produces a relatively hard spectrum, and for the value of oe and 
the cuts that we use, SPS tends to dominate over the full range of transverse momentum considered. On the other 
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Fig. 41. Event rate as a function of the angle between the two planes defined by the bb and jj systems. In SPS events, there 
is a correlation among the planes which is absent for DPS events. 


hand, DPS produces a much softer spectrum which (up to issues of normalization in the form of ogg) can dominate 
at small values of transverse momentum. The cross-over between the two contributions to the total event rate is ~ 30 
GeV for the acceptance cuts considered. A smaller (larger) value of cep would move the cross-over to a larger (smaller) 
value of the transverse momentum of the leading jet. 
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Fig. 42. The transverse momentum pr distribution of the leading jet in jjbb after minimal cuts. 


We turn next to the search for variables that may allow for a clear separation of the DPS and SPS contributions. 
Since the topology of the DPS events includes two 2 > 2 hard scattering events, the two pairs of jet objects are 
roughly back-to-back. We expect the azimuthal angle between the pairs of jets corresponding to each hard scattering 
event to be strongly peaked near A¢;; ~ App ~ T. Real radiation of an additional jet, where the extra jet is missed 
because it fails the threshold or acceptance cuts, allows smaller values of A@;;. The relevant distribution is shown for 
light jets (non b-tagged) in Fig. (43h. There is a clear peak near A¢;; = m for DPS events, while the events are more 


50 P. Bartalini et al.: Multi-Parton Interactions at the LHC 


broadly distributed in SPS events. The secondary peak near small Ad@;; arises from gluon splitting which typically 
produces nearly collinear jets. The suppression at still lower A¢;; comes from the isolation cut AR,; > 0.4. 
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Fig. 43. (a) The difference A¢ in the azimuthal angles of light jet pairs for DPS and both SPS+DPS events. The dijet pairs are 
back-to-back in DPS events. (b) The variable Sy for DPS and SPS+DPS events provides a stronger separation of the underlying 
DPS events from the total sample when compared to A¢ for any pair. 


_ The separation of DPS events from SPS events becomes more pronounced if information is used from both the 
bb and jj systems. As an example, we consider the distribution built from a combination of the azimuthal angle 
separations of both jj and bb pairs, using a variable adopted from Ref. [4]: 


pa JVA) FT. (67) 


In Fig. |43b, we present a distribution in Sg for both DPS and SPS+DPS events. Again, as in the case of the Ag 
distribution, the SPS events are broadly distributed across the allowed range of Sg. However, the combined information 
from both the bb and jj systems shows that the DPS events produce a sharp and substantial peak near Se ~ m which 
is well-separated from the total sample. 

The narrow peaks near A¢;; = m in Fig. |43p and near Sy = 1 in Fig. will be smeared somewhat once soft 
QCD radiation and other higher-order terms are included in the calculation. 

Another possibility for discerning DPS is the use of the total transverse momentum of both the bb and jj systems. 
At lowest order for a 2 —> 2 process, the vector sum of the transverse momenta of the final state pair vanishes. In 
reality, radiation and momentum mismeasurement smear the expected peak near zero. Nevertheless, the DPS events 
are expected to show a reasonably well-balanced distribution in the transverse momenta of the jet pairs. To encapsulate 
this expectation for both light jet pairs and b-tagged pairs, we use the variable [4]: 


lpr (bi, ba) ) | ( Ipr (1s j2) | (68) 
Sox -7 (mm (b1)| + ler)  \ pri) + lpr Ge)” 
Here pr(bj, b2) is the vector sum of the transverse momenta of the two final state b jets, and pr(j1, j2) is the vector 
sum of the transverse momenta of the two (non b) jets. 

The distribution in Dae is shown in Fig. |44| As expected, the DPS events are peaked near Spy ~ 0 and are well- 
separated from the total sample. The SPS events, on the other hand, tend to be far from a back-to-back configuration 
and, in fact, are peaked near Spr ~ 1. This behavior of the SPS events is presumably related to the fact that a large 
number of the bb or jj pairs arise from gluon splitting which yields a large pr imbalance and, thus, larger values of 
Spr 

Our simulations suggest that the variable S,„ may be a more effective discriminator than Sy. However, given the 
leading order nature of our calculation and the absence of smearing associated with initial state soft radiation, this 
picture may change and a variable such as Sy (or some other variable) may become a clearer signal of DPS at the 
LHC. Realistically, it would be valuable to study both distributions once LHC data are available in order to determine 
which is more instructive. 


P. Bartalini et al.: Multi-Parton Interactions at the LHC 51 


7000 
8000 


6000 


SPS+DPS SPS+DPS 


6000 5000 


4000 


4000 
3000 


N_events / Bin 


N_events / Bin 


2000 


oles a ek aF e E a ake 


Fig. 44. Distribution of events in Soe for the DPS and SPS samples. Due to the back-to-back nature of the 2 — 2 events in 
DPS scattering, the transverse momenta of the jet pair and of the b-tagged jet pair are small, resulting in a small value of Sp,- 
In (a) we show the Sp, distribution for our standard cuts, and in (b) we increase the cut on the transverse momentum of the 
leading jet, pi > 40 GeV. The fraction of DPS events in the whole sample decreases with increasing pr. 


The evidence in one-dimensional distributions for distinct regions of DPS dominance prompts the search for greater 
discrimination in a plane represented by a two dimensional distribution of one variable against another. One scatter 
plot with interesting features is displayed in Fig. ie DPS events are seen to be clustered near Sj, = 0 and 
are uniformly distributed in . The SPS events peak toward Spr = 1 and show a roughly sin® character. While 
already evident in one-dimensional distributions, these two features are more apparent in the scatter plot Fig. 
Moreover, the scatter plot shows a valley of relatively low density between Son ~ 0.1 and ~ 0.4. In an experimental 
one-dimensional @ distribution, one would see the sum of the DPS and SPS contributions. If structure is seen in data 
similar to that shown in the scatter plot Fig. |45| one could make a cut at Sp, < 0.1 or 0.2 and verify whether the 
experimental distribution in @ is flat as expected for DPS events. 


5.1.2 Strategy and Further Work 


The clear separation of DPS from SPS events in Fig. [45] suggests a methodology for the study of DPS. One can begin 
with a clean process such as pp —> 6bj1j2X and examine the distribution of events in the plane defined by ae and 
P. We expect to see a concentration of events near Sor = 0 that is uniformly distributed in &. These are the DPS 
events. Assuming that a valley of low density is observed between a ~ 0.1 and ~ 0.4, one can make a cut there that 
produces an enhanced DPS sample. Relative to the overall sample, this enhanced sample should show a more rapid 
decrease of the cross section as a function of the transverse momentum of the leading jet, and the enhanced sample 
can be used to measure Cep. A similar examination of other final states, such as 4 jet production, will answer whether 
the extracted values of oe are roughly the same. Theoretical and experimental studies of other processes can follow, 
such as bbtt, Wjj, and Hj. 

On the phenomenological front, next-to-leading order (NLO) expressions should be included for both the SPS and 
DPS contributions. The NLO effects are expected to change normalizations and, more importantly, the distributions 
in phase space. The sharp peaks near Ad;; = 7 in Fig. pah, Sg ~ 7 in Fig. [43p, and Spr =0Q0in Fig. [44] will be broader 
and likely displaced somewhat. 

Finally, it would be good to examine the theoretical underpinnings of Eq. and, in the process, gain better 
insight into the significance of ceg. A firm basis is desirable for Eq. (66) starting from the formal expression for the 
differential cross section in terms of the absolute square of the full matrix element integrated over phase space: 


>. 1 =.. 
do(pp — bbj1j2X) = 55 MP > bbjij2 X)| dP Sj jx: (69) 


The amplitude M (pp > bbjıj2 X) should include a sum of amplitudes for 2-parton collisions (one active from each 
incident hadron, i.e., 2 — 4); 3-parton collisions (two active from one hadron and one active from the other); and 
4-parton collisions (two active from each hadron or three from one and one from the other), and so forth that all yield 
the same 4 parton final state. There will be contributions to the final state from the squares of individual amplitudes 
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Fig. 45. Two-dimensional distribution of events in the inter-plane angle @ and the scaled transverse momentum variable Sp, 
for the DPS and SPS samples. 


as well as interference terms. Specializing to 4 + 4, the DPS case, one would start from a 4-parton — 4-parton hard 
part. Not evident at this time is how the four-parton matrix element can be reduced to a product of two matrix 
elements for the single parton scatterings, needed for Eq. (66). The demonstration of clear DPS signals in LHC data 
would be an important stimulus for further theoretical studies. 


5.1.3 NLO DPS Study of pp + WbbX — lvbbX 


In a subsequent paper [158], we investigate the possibility to observe double parton scattering at the early LHC in 
the process pp > WbbX — fvbbX. Our analysis is done at next-to-leading order in QCD. It begins with the basic 
assumption that Wbb production consists of two components: the traditional single parton scattering (SPS) process 
and the double parton scattering (DPS) process where two individual hard scatterings produce the Wb final state. 

After identifying the most relevant background processes, we pinpoint a set of observables and cuts which would 
allow for the best separation between the DPS Wb signal and the backgrounds (including the SPS Wbb process). To 
provide the most precise predictions possible, we generate the DPS Wbb signal event sample, the SPS Wbb sample, 
and the dominant background event samples at next-to-leading order in QCD. The main obstacles in the extraction 
of the DPS signal are the backgrounds from tt production and the SPS Wbb component. The most efficient way to 
suppress the tt background is with an upper cut on the missing transverse energy of the event, since top quark decays 
result in larger values of Fr. 

To separate the DPS component of Wbb from the SPS component, we find it useful again to employ observables 
which take into account information on the full final state rather than observables which involve one or two particles. 
Examples are the Sj, variable and the angle (AO,; »,) between the two planes defined by the bb and £v systems, 
respectively. By displaying the information from these two observables in two-dimensional distributions, we show that 
it is possible to identify distinct regions in phase space where the DPS events reside. Utilizing cuts on these observables 
that enhance the DPS Wbb sample, we find that the DPS signal can be observed with a statistical significance in the 
range S/V B ~ 12 — 15. 

The focus in Ref. [158] is to establish double parton scattering as a discernible physics process at LHC energies 
and measuring the size of its contribution. Once DPS production of Wbb is observed, it will be interesting to assess its 


P. Bartalini et al.: Multi-Parton Interactions at the LHC 53 


TGS09 OMSTWo OMSTW, _OMSTW2 
W'W | 0.546 0.496 0.409 0.348 
Wtw | 0.321 0.338 0.269 0.223 
WW | 0.182 0.182 0.156 0.136 
R 
0.784 1.00 1.00 1.00 


Table 2. DPS WW total cross sections (in pb) and ratio R defined in Eq. for pp collisions at ys = 14 TeV evaluated 
using different dPDFs sets. 


OGso9 
vs =7 TeV vys= 10 TeV vys= 14 TeV 
wtw 0.107 0.250 0.546 
wrwr 0.0640 0.148 0.321 
Ww-w- 0.0317 0.0793 0.182 
R 
0.709 0.751 0.784 


Table 3. DPS WW total cross sections (in pb) and ratio R defined in Eq. (70) for pp collisions at different CM energies y's. 


potential significance as a background in searches for other physics, such as Higgs boson production in association with 
a W boson (where the Higgs boson decays as H — bb), and precise studies of single top quark production where new 
physics could contribute to the Wtb vertex. A detailed analysis of either of these channels would require a different 
set of optimized physics cuts. We limit ourselves in Ref. to showing the bb invariant mass distribution for the 
évbb final state. We see that the DPS Wbb component alters the overall shape of the bb mass spectrum, enhancing the 
small mass region. This feature is consistent with our earlier observation that the pr spectrum of leading jets is softer 
in the DPS component. The DPS component contributes primarily in the region below 120 GeV or so. At face value, 
it does not seem to pose a hindrance for searches for Higgs bosons in the HW channel. However, Wbb DPS could 
be a significant background in the search for new particles, with masses in the 50 - 100 GeV range and appearing as 
resonances in M, and it should be accounted for in any analysis. 


5.2 Probing double parton scattering with same-sign W pairs at the LHC 
Contributing authors: J. R. Gaunt, C.-H. Kom, A. Kulesza, and W. J. Stirling 


In this section we focus on W+W= production, followed by decays into same-sign di-leptons (SSDL) plus missing 
transverse energies (Ær). Our discussion is based on (see for an early study of this process). One advantage 
in studying DPS W+W-= is that single WF production is well described theoretically, and will be accurately measured 
at the LHC. As a result an approximation with no correlations included will be well modelled. In addition, W+W 
production through single parton scattering (SPS) is forbidden at the same order in the SM, i.e. there is no gq > 
W+W= contribution, due to (U(1)gm) charge conservation. The lowest order ‘background’ process is instead qq’ > 
W+W+*jj, which is of order O(a*) or O(a2,a7). The additional jet activity can be used to distinguish this background 
from the DPS signal. 

Same-sign W*+W-= production could then potentially be a benchmark DPS process. We will study kinematic 
properties of this process by including longitudinal correlation effects, namely constraints from momentum and valence 
number conservation. These constraints are consistently included in the double parton distribution function (dPDF) 
set GS09 [82]. We refer readers to Section [4.3] and original literature for details on the construction of this dPDF set. 
We will quantitatively study effects of these constraints by comparing results using GS09 with other simple dPDF 
models that will be defined later. In the literature, comparison of DPS and SPS W*+W*= production has mostly been 
studied at the level of production cross sections. Here we discuss strategies to extract the signal from the background 
including important contributions from heavy flavour and gauge boson pair productions. 


5.2.1 The W+W*= signal 


Our phenomenological investigation for DPS W*+W*= production is based on the equation [64] with m = 1, A= B = 
W= and tı = t2 = ln( Mẹ). We assume ceg = 14.5mb, the value obtained by the CDF y + 37 B] analysis. We adopt 
GS09, which incorporates momentum and sum rule constraints, as the default longitudinal dPDF choice. To see the 
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Fig. 46. Normalised /* pseudorapidity distributions for pp collisions at ys = 14 TeV evaluated using different dPDFs. No cuts 
are applied. 


effect of these constraints, we compare cross sections using factorised dPDFs of the form: 


Djj(@1, £2) = Di(£1)D; (x2) 
xO(1— zı — £2)(1 — z1 — 22)", 
m= 0,1,2, (70) 


hereafter referred to as MSTW, sets. In the above expression, the factorization scale is again fixed at up = Mw. The 
completely factorised approximation that has been taken in existing phenomenological studies is obtained by setting 
D;j(£1, £2) = D;(x1)D;(x2). The sPDFs D;(x) are taken from the MSTW 2008 LO set. In Table [2] we compare the 
cross sections and the values of the ratio 


Ras WOWW, (71) 
Ow+w- 


which measures the deviation from the factorisation approach (R = 1 when factorisation is exact) using different 
dPDFs. In Table [3| the value R as a function of CM energy using GS09 is evaluated. We see that using GS09, 
factorisation is broken at the 20% to 30% level, and the approximation improves at the higher collider energies as 
lower x regions are probed. On the contrary, the MSTW,, models have R values very close to 1 because effects of the 
momentum suppression factors (1 — zı — x2)” in the numerator and the denominator of the expression tend to cancel. 

Sum rule effects incorporated in GSO9 can also be seen in the distributions of the charged leptons [+ from W= 
decays. To understand this, note that production of same sign W~’s along the same forward direction is suppressed 
compared with factorised models due to both momentum and number sum rule constraints. This is because this phase 
space region is significantly influenced by D”” (x1, £2), where v = vu, va denotes valence up or down quarks, both of 
which with relatively high momentum fractions x; and x2. In fact, D”4”4 = 0 identically in GS09 due to the presence 
of only one valence down quark in a proton, whereas it is non-zero for other simple factorised models. The /* hence 
have more central distributions compared with the MSTW,, models. The pseudorapidity (m) distributions for I+ are 
shown in Fig. 

Such longitudinal correlations also imply that the SSDL pairs prefer to lie in opposite hemispheres (i.e. m, Xm, < 0). 
This preference can be quantified by computing pseudorapidity asymmetry, defined as 


= o(m, X Ms < 0) — a(n, X Me > 0) 


a 7 
i a(n, X Na < 0) T a(n, X My > 0) 


(72) 


where |m, |, |n| > 7". The value of ap for [+ as a function of n™™ is displayed in Fig. It increases with 7j™™ 
as the correlations are most important for the distributions probed at high values of x for both partons in the same 
proton, which is reached when the leptons are produced at high |r|. For GS09, this effect is even more pronounced for 
I~, as only one valence down quark is present in a proton, making simultaneous extraction of two high x down quarks 
highly suppressed. 
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Fig. 47. Pseudorapidity asymmetry a,, for pp collisions at ys = 14 TeV evaluated using different dPDFs. No cuts are applied. 
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Fig. 48. Cross sections of various electroweak processes in p-p(p) collisions as a function of ys. The dotted curves correspond 
to single scattering processes, while the solid curves correspond to double scattering processes computed using GS09 dPDFs. 


5.2.2 Single scattering backgrounds 


We now turn to the single scattering backgrounds. We show in Fig. [48] the total cross sections of the DPS signal, as well 
as electroweak single scattering backgrounds and related processes. As discussed before, due to charge conservation 
the lowest order ‘irreducible’ SPS background including a pair of same sign W’s is q? — W*W=jj. At the 14 TeV 
LHC, the DPS signal and SPS background total cross sections are of the same order of magnitude (O(0.1—1)pb before 
leptonic branching ratios). However the presence of two additional hard jets means that central hard pr jet veto can 
be useful in suppressing this background. The SPS cross section as a function of max jet pr and min jet 7 at the 
parton level is displayed in Fig. As we can see, this background can be suppressed easily, for instance by vetoing 
events with central jets with pr > 20 GeV. 
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Fig. 49. o(pp > WTW *jj)- [BR(WF > ptv,)]? (fb) as a function of max pr; and n?™™. No other cuts are applied. 
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Actually other electroweak processes that can lead to the same final states. In particular, both W*Z(7*) and 
Z(y*)Z(y*) gauge boson production can lead to SSDL + Ær when the ‘wrong’ sign leptons are not identified, for 
instance if they lie outside the central tracking region, or if they are too soft to be reconstructed. The relevant 
processes are 


qq > WIZ(y*) > vl (Ll), 
qg > Z(O) ZO“) oO), (73) 


and their charge-conjugated processes. In the above expressions, charged leptons in brackets are not identified. 

As can be seen in Fig. the total cross sections for W*Z and ZZ are about an order of magnitude above the 
signal, while the cross sections are significantly larger when y* is involved. However, for these processes the lepton pr 
spectra are much harder compared with the signal, hence a max lepton pr cut can reduce this background. A wrong 
sign lepton-veto in the central region will also be useful in reducing the Z contribution, while looking for presence of 
low invariant mass system of an isolated charged tracks and a nearby identified lepton [160] can help suppressing the 
7* contribution. 7 

Another source of background is pair production of heavy flavour quarks. Production of bb can lead to SSDL pairs, 
when a neutral B-meson is present and undergos B°-B° mixing, followed by semi-leptonic decay for both B-mesons. 
The relevant processes are then: 


gg > bb > BB+..., 
B > ItvX, 
B° > B° yx, (74) 


together with the charge conjugation processes. 

The bb cross section is orders of magnitude larger than that of the signal. However, the pr spectrum of the B-mesons 
decreases exponentially. It is thus very difficult for both charged lepton and neutrino from the semi-leptonic B decay 
to acquire large transverse momenta. Imposing tight lepton isolation, min lepton pr and min Ær cuts will thus help 
suppressing contribution from this process. 

Production of tt pair can also result in SSDL pairs when a top and the bottom of the other top decay semi- 
leptonically. The relevant processes are 


t => Wtb > l*vb, 
t> Wb > qg'l*ve. (75) 


These events have significant jet activities. Also, the lepton pr spectrum is much harder compared with the signal. 
The leptons, particularly the ones from B-decays, are usually poorly isolated. We thus expect tight lepton isolation, 
central jet veto and a max lepton pr cut will be effective in suppressing this background, and will not consider this 
process further in our numerical analysis. 
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| Outpt (fb) Ou pe (fb) 
| W*W~(DPS) 0.82 0.46 
W*Z(9*) 5.1 3.6 
Z(¥"*)Z(7*) 0.84 0.67 
bb (py > 20 GeV) 0.43 0.43 


Table 4. Cross sections (in fb) of the processes simulated after cuts, including branching ratios corresponding to same-sign 
dimuon production. 


5.2.3 Observing DPS W+W*= at the LHC 


We see in the last section, that despite the backgrounds having much larger cross sections than the signal, there are 
handles to suppress the former processes. To see to what extent the signal can be extracted from these backgrounds, 
we perform a parton level signal + background simulatior}| with the following criteria: 


— Both leptons in the like sign lepton pair must have pseudorapidity |ņ| < 2.5. 

— Both leptons are required to be isolated: 

El o < Ez™ = 10 GeV, where El o is the hadronic transverse energy in a cone of R = 0.4 surrounding each of the 
like-sign leptons. 

— The transverse momenta of both leptons, p/,, must satisfy 20 < pl, < 60 GeV. 

— An event is rejected whenever a third, opposite-signed, lepton is identified. A lepton is assumed to be identified 
with 100% efficiency when ph > pid and |n| < n't, where pid = 10 GeV and n" = 2.5. 

— The missing transverse energy Ær of an event must satisfy Zr > 20 GeV. 

— Reject an event if a charged (lepton) track with pis > pr > 1 GeV forms an invariant mass < 1 GeV with one of 
the same-sign leptons. 


The cross sections after these cuts are displayed in Table[4| The largest background comes from the W* Z(y*) pro- 
cesses, which is a factor of a few larger than the signal. Unfortunately, we find that many basic kinematic distributions 
are fairly similar between this SPS background and the DPS signal, making further cuts unlikely to be beneficial. 

On the positive side, it might be advantageous to exploit the fact that the value of ay, is relatively small, but positive, 
for the DPS signal. On the contrary, a,, tends to be negative for the background, which reflects the preference for the 
leptons to lie close in pseudorapidity space in order to reduce the CM energy of the system. This property is illustrated 
in the diagram on the left of Fig. which shows how a», for different processes vary as a function of a minimum m 
cut. 

The ratio of positively charged (++) and negatively charged (——) SSDL events (which we call charge asymmetry 
ratio) may also be used, as initial state partons of various flavours and momentum fractions are involved in different 
processes, each leading to different charge asymmetry ratio. The charge asymmetry ratio for the DPS signal and 
SPS W+Z(y*) background using different lepton identification criteria is displayed in the right figure of Fig. An 
important point that should be noted is that this ratio is fairly stable when varying the cuts, as can be inferred from 
the results in Table [4] 

To summarise, same-sign W*W-~ production can potentially be a benchmark process for studying double parton 
scattering at the LHC. We have discussed likely changes to the kinematic properties of the signal when momentum 
and sum rule constraints are included in the description of dPDFs. We find that, after including important physics 
backgrounds previously overlooked in the literature, a small excess of SSDL DPS events could be observed at the LHC, 
while further improvements can be made to enhance the signal. Finally, we note that there are other DPS processes 
with purely leptonic final states, such as double Drell-Yan and double J/psi production, that could be interesting. 
Recent phenomenological studies of these processes can be found in Refs [161|[162)(763)/164)/165) . 


5.3 Multiple Parton Interactions in Z + jets production at the LHC 
Contributing author: E. Maina 


5.3.1 Multiple Parton Interactions in Z + jets production at the LHC 


In Ref. [166] the contribution of MPI to Z + n-jets production at the LHC, n = 2,3,4, where the Z boson is assumed 
to decay leptonically, has been examined. These processes have the advantage of a much larger cross section than 
same-sign WW production and therefore are more likely to allow detailed studies of MPI at the low luminosity, about 
1 fb~', foreseen for the first two years of operation at the LHC with ys = 7 TeV. 


3 We refer readers to the original paper for a technical account of the simulation. 
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Fig. 50. Left: pseudorapidity asymmetry a,, for the positive SSDL+ẸÆr DPS signal and selected SPS background, after imposing 
cuts described in the text. Right: charge asymmetry ratio (++)/(——) as a function of lepton identification criteria for different 
processes. 


Z + nj production probes initial state parton combinations which are different from those probed in WSW 
processes. The latter, at lowest order, are always initiated by four-fermion states, mainly udud. The former, on the 
contrary, typically have at least two gluons in the initial state since the largest component [167168] involves a two jet 
process which is dominated by gluon-gluon scattering. 

For comparison we also present the predictions for y + 3j production, the reaction from which the most recent 
and precise estimates of o,f have been extracted at the Tevatron. The CDF collaboration [23] has measured eff = 
14.5 + 1.7154 mb, a value confirmed by DO which quotes coeff = 15.1 + 1.9 mb Ø]. In Ref. [79] it is argued, on the 
basis of the simplest two channel eikonal model for the proton—proton cross section, that a more appropriate value at 
vs = 1.8 TeV is 10 mb which translates at the LHC into one = 12 mb. Treleani then estimates the effect of the 
removal by CDF of triple parton interactions (TPI) events from their sample and concludes that the CDF measurement 
yields o-¢ ~ 11 mb at Tevatron energies. In the following we use o-+¢¢ = 12.0 mb for all LHC center of mass energies, 
with the understanding that this value is affected by an experimental uncertainty of about 15% and that it agrees 
only within 30% with the predictions of the eikonal model. Since geff appears as an overall factor in our results it is 
easy to take into account a different value. 

It is worth mentioning that at present there is a discrepancy between the value of oes extracted by CDF and 
DO and the one which is effectively employed by PYTHIA whose normalization is derived mainly from comparisons 
with small pr data which dominate the total cross section. The description of MPI in PYTHIA 8 [56] assumes that 
interactions can occur at different pr values independently of each other inside inelastic non—diffractive events. The 
expression for a DPS cross section, here referred to as Double Parton Interactions (DPI) cross section, becomes 
therefore: 


o=< fimpa > O° o2/onp/k (76) 


where oyp is the total non-diffractive cross section and fimpact is an enhancement/depletion factor chosen event- 
by-event to account for correlations introduced by the centrality of the collision. This quantity is typically averaged 
during an entire run to calculate < fimpact > in Eq. Typical values at the center of mass energy of 10 TeV are 
1.33 for < fimpact > and 51.6 mb for oy p. Comparing Eq. with Eq. [65] tells us that PYTHIA 8 predicts an effective 
Sef f=InD/< fimpact > Which is about a factor three larger than the one actually measured at the Tevatron. I believe 
that this issue deserves careful consideration and that new measurements of high pr MPI reactions would be quite 
welcome. 

NLO QCD corrections are or will soon be available for all SPS, here referred to as Single Parton Interactions (SPI) 
processes leading to an electroweak vector boson in association with up to four jets [L69)(T70\1 711 72)(1 73/1 74/1 7511 76) . 
The Drell-Yan cross section is known at NNLO [177]. Measurements at the Tevatron show good agreement between 
NLO calculations and data [T78|[179]. These new developments open the possibility of validating the predictions using 
events with large visible energy, where the MPI contribution is small, and then using them for a direct measurement 
of the MPI cross section at smaller total invariant masses in parallel with more data driven analysis similar to those 
of CDF and DO. 

In the following we compare the results obtained with the GS09 dPDF with those obtained with two instances 
of fully factorized sPDF: MSTW2008LO and CTEQ6L1 [40]. Hence we can estimate, even in the absence of a 
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proper dPDF set based on CTEQ6, the dependence of MPI predictions on the choice of PDF, a study that to our 
knowledge has not been performed before. 
The MPI processes which contribute at leading order to Z + n—jets through Double Parton Interactions are those 
in which an event producing k jets is superimposed to an event producing a Z-boson and (n — k) jets, k = 2,...,n. 
All samples have been generated with the following cuts: 


pr, > 30 GeV, |m] <5.0, 
pr, = 20 GeV, |ne] < 2.5, (77) 
pr, = 30 GeV, |n| < 2.5, 

AR; >0.5, ARy>01, AR; > 01 


where j = u,u,d,d,s,5,c,é,b,b,g andl = et,e7, u7, pr. 

The Z+4-jets sample has been generated with PHANTOM [ISOISIII82], while all other samples have been produced 
with MADVENT [1831184]. All samples have been generated using CTEQ6L1 parton distribution functions. The 
QCD scale (both in a, and in the parton distribution functions) has been taken as 


P= phi, (78) 
{=I 


where n is the number of final state partons, for all reactions with the exception of gg > [17 for which the scale has 
been set at Q? = M3. 

The results shown in the following under the CTEQ heading have been obtained combining at random one event 
from each of the reactions which together produce the desired final state through MPI. 

The results shown under the MSTW and GS09 headings have been obtained through a reweighting procedure 
by the appropriate ratio of parton distribution functions and coupling constants. For instance, an event like (qq; > 
gl*l-) ® (gg > gg), constructed from two events generated separately with CTEQ6 PDF, can be transformed in a 
weighted event with MSTW2008 PDF multiplying its original weight by 

go Ds” (tı) a" (tx) p" wm (t2) Fa (t2)? 


R=— 


x x x 
DEF? (DT) OST (th) ~ DOF (ta) DS?) a (to)? 


(79) 


where tı, t2 are the factorization scales for qg@ —> gltl~ and gg — gg respectively. The factorization scales have 
been read off from the event files. The second and fourth factors in Eq.(79) take into account the different values of the 
strong coupling constants for the two different sets of PDF: äro (Mz) = 0.130 while asto (Mz) = a’ (Mz) = 
0.139. The only difference for the GS09 case would be that the correlated dPDF Fi; (t1,t2) would appear instead of 
the uncorrelated product D;(t:)D,(t2) and so on. 

All results are obtained with the following values for the electroweak input parameters: Mz = 91.188 GeV, My 
= 80.40 GeV, Gr = 0.116639 x 1075 GeV~?. 

The total cross sections for SPI and DPI production for Z + 2-jets, Z + 3-jets and Z + 4-jets are presented in 
Table [5] In our estimates below we have only taken into account the muon decay of the Z boson. 


14 TeV 10 TeV 7 TeV 
CTEQ | MSTW | GS09 | CTEQ | MSTW | GS09 | CTEQ | MSTW | GS09 
Z+2) 
SPI | 52.65 | 60.70 30.63 | 35.15 16.56 | 18.88 
DPI | 11.27 | 1437 | 15.50 | 4.80 6.35 | 6.68 | 1.88 261 | 2.66 
Z +3) 
SPI | 15.71 | 19.10 8.46 | 10.23 411 4.93 
DPI | 2.70 3.75 | 3.88 | 1.02 149 | 148 | 0.34 0.54 | O51 
Z+4j 
SPI | 4.26 5.41 2.00 2.53 0.83 1.04 
DPI | 0.96 1.53 | 150 | 0.33 0.56 | 0.52 | 0.10 0.18 | 0.16 


Table 5. Z + n-jets, Z > u™ u” cross sections in pb. Cuts as in Eq. (77) with AR;; = 0.5. 


The total cross sections for SPI and DPI production for y + 3-jets are shown in Table |6} At the LHC trigger 
thresholds for single photons are foreseen to be much higher than those for double leptons |186//187\|188}. While pair 
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of leptons are expected to be triggered on for transverse momenta of about 15 GeV, single photons will be detected 
only when their transverse momenta is larger than about 80 GeV at the design energy of 14 TeV. Since MPI processes 
are known to decrease sharply with increasing transverse momenta, we present also the predictions for pr, > 80 GeV. 

The Single Particle Interaction MSTW results are larger than those obtained with the CTEQ PDF by an amount 
which varies between 15% for Z +27 to 27% for Z + 4j, increasing as expected with the power of a, in the amplitude. 
The Double Particle Interaction MSTW results are larger than those obtained with the CTEQ PDF by an amount 
which varies between 30% and 90%. The larger shift is due to the smaller scales for the two individual scatterings 
compared to a single interaction event with the same final state particles. The predictions for the GS09 correlated 
dPDF are larger than those with MSTW uncorrelated ones for ys = 14 TeV and ys = 10 TeV while they are smaller 
for ys = 7 TeV. The difference is at most of 15%. Taking into account the errors in the measurement of off we 
conclude that the uncertainties due to the choice of PDF and to correlation effects are reasonably under control. 

These variations should be compared with the uncertainty due to scale variation in PDF and in the strong coupling 
constant. Changing the scale in Eq. (78) by a factor of two in either direction for two limiting cases, namely Z + 27 
production at ys = 7 TeV and Z + 47 production at \/s = 14 TeV the cross section changes by +14%/-13% in the 
first case and by +57%/-29% in the second. 

The effects of higher order corrections are more difficult to estimate since no NLO calculation for MPI processes is 
available. QCD one loop calculations are available for vector boson production with up to four jets and are typically of 
order 10% with the exception of Drell-Yan inclusive production where they are of the order of 50%. NLO corrections 
for the inclusive jet cross section at the LHC have been presented in Ref. [189]. For small transverse momenta, as the 
ones we are interested in this paper, they are of the order of 10%. 

The ratio between the MPI and SPI cross sections increases with the collider energy, that is with decreasing average 
momentum fractions carried by the incoming partons. It also increases with the AR;; separation because of the absence 
of correlations between the final state partons originating in the independent scatterings which compose MPI events. 
For Z+nj processes and taking AR,;; = 0.5 as an example, the ratio is of the order of 10% for ys = 7 TeV and grows 
to about 25% at \/s = 14 TeV. The results for y + 3-jets show a similar behaviour with somewhat smaller fractions 
of MPI events to SPI ones which however depend drastically on the pr, cut. For pr, > 30 GeV they range between 5 
and 10% while for pr, > 80 GeV they are at the percent level. 


14 TeV 10 TeV 7 TeV 
y+3j | CTEQ | MSTW | GSO9 | CTEQ | MSTW | GS09 | CTEQ | MSTW | GS09 
pr, = 30 GeV 
SPI 4516.7 | 5610.2 2637.2 | 3263.8 1415.8 | 1744.6 
DPI 422.2 593.7 642.2 | 170.3 254.1 264.9 62.0 100.0 99.4 
pr, = 80 GeV 
SPI 671.5 813.0 368.30 | 443.38 177.4 212.0 
DPI 17.7 28.1 28.1 6.59 11.49 10.85 2.09 4.09 3.58 
Table 6. y + 3-jets cross sections in pb. Cuts as in Eq.{77p. 


If we consider the MPI processes as our signal and the SPI ones as the corresponding background, we can estimate 
the prospect of measuring MPI in a given final state from the standard S/V B significance. Using for S the result 
obtained with GS09 PDF and for B the result for the MSTW set and assuming a luminosity of one inverse femtobarn 
at 7 TeV, the significancies extracted from Table [5] in the Z + ptp- channel alone, are 19/7/5 for 7+2/3/4 jets . 
The corresponding number of expected MPI events are 2600/500/160. Therefore it appears quite feasible to measure 
the MPI contribution to 7+2/3/4 jets already in the first phase of the LHC. 

The significance of y + 3-jets depends on the trigger strategies. If the threshold for single photon detection can 
be brought in the 30 GeV range then the much larger production rate, about ten times that of Z(uj4) + 27, provides 
the best opportunity for an early measurement of MPI at the LHC. If, on the contrary, the photon trigger cannot 
substantially deviate from about 80 GeV, Z + 2 production looks more promising than the y + 3—jets channel whose 
significance becomes similar to that of Z + 37. 

The contribution to the MPI Z + n-jets cross section due to two jet production in association to Z + (n — 2)-jets 
processes is in all instances the largest one, therefore, even with more than two jets in the final state, the majority of 
MPI events are expected to contain a pair of jets which are back to back in the transverse plane. This is confirmed by 
the left hand side of ea which displays the distribution of the angular separation A¢ between the two highest pr 
jets in Z + 4j events at ys = 7 TeV and AR;; = 0.7. 

The right hand side of Fig. presents the total visible mass distribution in Z + 27 production with the same 
energy and angular separation. It clearly shows that MPI events are produced with a smaller center of mass energy 
than SPI ones. 
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Fig. 51. On the left: distribution of the angular separation in the transverse plane between the two highest pr jets in Z + 47 
events. On the right: distribution of the total visible mass, ()7"_, pi)”, in Z +24 events. For both plots ys = 7 TeV, ARj; = 0.7. 


6 Summary and Outlook 


The understanding multi-parton interactions (MPI) in hadronic collisions remains a challenge. Traditionally, research 
has concentrated on four tasks covered during this workshop: experimental measurements of underlying events and 
minimum bias events, models of soft physics implemented in Monte Carlo (MC) generators, development of the 
theoretical description of MPI, and phenomenological studies. In particular, a lot of effort, strongly driven by the 
experimental measurements, reported here by the ATLAS, CMS and LHCb collaborations, see Sections [2.1] 
respectively, focuses on modeling of minimum bias and underlying event physics in MC generators. This is, of course, 
only natural and very much needed given that the MC generators are major tools used in the experimental analyses 
by the LHC collaborations. In fact, one of the main tasks in LHC phenomenology is the construction and development 
of general purpose MC generator fully able to describe exclusive states, including contributions from MPI. This task 
has enormous importance for the overall success of the LHC physics programme. However, in order to improve the 
simulation of the MPI in the MC codes, more elaborate theoretical input, supplemented by information from high-pr 
phenomenological studies, is necessary. 

The standard MC algorithms are based on the factorization theorem and the factorization of the QCD amplitudes. 
This theorem lets us to define general algorithms and leads to good predictions for the LHC. However, the extension 
of the single interaction picture to encompass MPI is not a simple task. The main challenge is to understand the 
structure and the topology of the MPI based on perturbative QCD, as discussed by D. Treleani in Section [4.4] Two 
important issues can be identified: the treatment of the parton distributions and the parton dynamics. 

The presence of multiple interacting partons in the initial state enforces the appearance multi parton distribution 
functions (mPDF). In the MC programs the mPDF functions are simple scalar functions, resulting from approximating 
a mPDF function with n partons by a product of single PDF functions. This is a good approximation under the 
assumption that the correlations in colour, spin and flavour can be neglected, as the most important information about 
the mPDF in the parton shower context comes from their evolution. Deficiencies in the approximation of the initial- 
state mPDF reached after backward evolution can be compensated through changes (tuning) in the hadronization. 

In the current MC implementations the evolution of the multiple partons taking part in the same collision relies 
on well understood evolution equations for single PDF. Going beyond this approximation and moving to full mPDF 
evolution poses many challenges. The evolution equation for mPDF scalar functions has been available for some time. 
Recently, J. R. Gaunt and W. J. Stirling, cf. Section provided a solution and detailed studies for the double 
parton case. However, the question of how double parton distributions should be defined and which evolution equation 
they satisfy has been reconsidered since the workshop, see [9590191]. Moreover, it is not clear whether the mPDF 
functions should be scalar functions or operators in color, spin and flavor space. As M. Diehl has pointed out in 
section [4.1] the non-trivial color, spin and flavor correlations occur in the dPDF case and they can be as important 
as the uncorrelated contributions. For example, non-trivial correlations in colour space can lead to similar effects as 
those caused by colour reconnection in the hadronization model. These issues would have to be resolved in order to 
achieve proper description of mPDF in MC generators. 

The other important ingredient in the MC models is the treatment of the parton dynamics. Initial and final state 
radiation is well understood in the standard parton shower models, but the MPI contributions are more complex. 
HERWIG++ (cf. Section follows a rather minimalistic approach and considers MPI effects only via the 2 > 2 
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parton interactions while PYTHIA studies also effects of rescattering and joint interaction in the MPI models, see 
Section These contributions are considered in a classical probabilistic framework. To better understand and 
improve the existing models it is important to look beyond the existing approach. In perturbative parton dynamics 
there are usually two types of contributions. The first type is real radiation, when one or more partons appear in 
the final state, and the partons are resolvable. These contributions can be described by a rather simple splitting 
function (operator) that can be obtained from tree level Feynman graphs. The other contributions are the virtual 
and unresolvable configurations which are not considered explicitly but are included via unitary conditions. Usually, 
the splitting operator of virtual contributions has some imaginary part that cannot be obtained from the unitary 
conditions and thus is not taken into account in the current MC implementations. Such imaginary contributions, 
known in the literature as the Coulomb gluon terms, can lead to the appearance of the so-called super leading 
logarithms in perturbative calculations for certain quantities. Correspondingly, these effects cannot be reproduced by 
the existing parton shower MC codes. Since the imaginary part of the virtual splitting operator can change the colour 
structure, its effect might be similar in the colour space to the action of colour reconnection. It has been pointed 
out that colour reconnection, contributing to hadronization, can have a big impact on the predictions, and therefore 
is important for description of data. Consequently, one should systematically consider all the effects of perturbative 
origin, e.g. the colour-changing effects, in order to ensure the universality of the non-perturbative hadronization model. 

It is often said that a good MC program has tunable parameters only in the hadronization model. This is more or 
less is true if we do not have to consider MPI effects, otherwise we need to deal with mPDF which have transverse 
momentum dependence. Unfortunately, we have very little information on the transverse momentum dependence. The 
standard practice is to assume the factorization of the longitudinal and transverse dependent parts of the mPDF. In 
the MC generators, the transverse momentum dependent part is described and parametrized by a simple function for 
every parton flavour. HERWIG++ makes use of only three tunable parameters while PYTHIA has more complicated 
parametrization with more parameters. This tunable function describes non-perturbative effects, and it can be shown 
that the average number of the interactions strongly depends on it. In principle we know very little about goodness 
of the factorization assumption. Theoretical studies, based on modeling the proton substructure by a dipole cascade 
model discussed by G. Gustafson in Section [4.5] show that the transverse part has a non-trivial dependence on the 
longitudinal momentum fraction of the incoming parton and the factorization scale. At this point, phenomenological 
studies of DPS processes with high pr final states have an important role to play, as comparisons with experimental 
results will serve to test the predictions and uncover further information on the transverse momentum dependence of 
the mPDF. 

The DPS measurements will also provide an important validation procedure of the MPI models built into MC gen- 
erators. For this purpose, one needs to identify the processes where the DPS signal is favourable over SPS background. 
Promising candidate processes such as same-sign W production, Z production in association with jets, four-jet pro- 
duction or i eee of a bb pair with two jets have been discussed in previous chapters. As pointed out by E. Maina 


in Section measuring DPS in different processes, for example in Z + jets and W*W~ production, will deliver 
complementary information since different initial state parton combinations will be probed. It is also important to 
study if one can define quantities more sensitive to DPS than the ones measured currently. Two examples of such 
studies for bb pair plus two jet production and pp ~ WbbX — lvbbX were presented by E. Berger, see Section|5.1| Al- 
ternatively, kinematical regions where the DPS provides the dominant signal can be selected. As reported by B. Blok, 
cf. Section in the case of four-jet production such region is constituted by back-to-back dijet production. For 
gauge boson pair production, the relevant region is where the transverse momentum of each boson is small, as shown 
in Section [4.1 It has to be also checked that the background to the actual final state observed experimentally can be 
sufficiently suppressed, cf. same-sign lepton final states discussed by C.-H. Kom in Section [5.2] 

In summary, understanding MPI in hadronic collisions requires further efforts on both theoretical and experimental 
fronts. In particular, it would be helpful to critically reevluate the description of MPI in MC codes and the tuning 
strategies. Although comparing results provided by existing MC codes supplemented with various tunes is certainly 
valuable on its own, there is a serious risk that it will not bring full understanding of the limitations and the sys- 
tematic errors of MC generators. Apart from purely theoretical work and the work related to development of MC 
codes, it would be advantageous to, for example, identify and explore phase-space regions sensitive to MPI effects, 
other than those used so far in the experimental studies. Furthermore, one would also benefit from constructing new 
variables probing MPI in particular processes. As already stressed, such efforts need to be undertaken in common by 
the experimental and theoretical communities. 
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